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Heat  transfer  models  are  effective  tools  for  the  analysis  of  thermal  processing  but 
their  application  is  generally  limited  to  simple  geometries  such  as  cylinders,  rectangles  and 
spheres.  A conical  and  a finite  elliptical  cylinder  are  alternative  shapes  for  products  to  be 
thermally  processed.  The  objectives  of  this  work  were  to  develop  heat  transfer  models  to 
predict  temperature  distribution  for  conduction  heated  foods  in  any  given  size  and  shape 
of  conically  and  elliptically  shaped  containers,  to  experimentally  validate  the  models,  and 


ity  or  nutrient  retention  under  different  processing 
conditions  with  the  goal  of  process  optimization. 

A finite  difference  method  with  energy  balance  approach  was  used  for  developing 
the  models.  The  conically  or  elliptically  shaped  containers  was  divided  into  small  volume 


elements.  Energy  balance  equations  were  developed  for  each  volume  element  and  solved 
explicitly. 

An  acrylic  material  of  known  thermal  properties  was  used  to  manufacture  different 
geometries  and  sizes  of  3 cones  and  2 finite  elliptical  cylinder.  Each  geometry  had  3 or  4 
thermocouples  (36  gauge,  T type)  inserted  at  different  locations.  Heat  penetration  tests 
were  carried  out  in  a water  bath  with  constant  or  variable  water  temperatures. 
Experimental  temperatures  agreed  well  with  temperatures  predicted  by  the  models. 

The  validated  models  were  used  to  calculate  lethality  and  quality  retention  under 
different  processing  conditions.  All  calculations  were  performed  for  a target  lethality  (F„) 
of  three  minutes  at  the  cold  point.  For  a conically  shaped  container,  the  location  of  the 
cold  point  was  searched.  The  effects  of  container  geometry , food  thermal  diffusivity 
(from  1 .2  to  1.6  mVs),  and  surface  heat  transfer  coefficient  on  the  retention  of  quality 
parameters  were  investigated.  The  kinetic  parameters  of  quality  attributes  were  the  D,i,°c 
value  (from  40  to  240  min)  and  the  Z value  (from  10  to  60  *C)  which  covered  the  range 
found  in  the  literature.  The  range  of  the  retort  temperature  was  from  1 10  to  14 OX  The 
foods  of  high  thermal  diffusivity  and/or  foods  packed  in  a container  leading  to  a rapid  heat 
transfer  to  the  cold  point  required  short  processing  times  and  retained  more  quality 
Within  the  scope  of  this  study,  the  optimum  temperature  for  maximizing  quality  retention 
was  highly  dependent  on  the  Z value  of  the  quality  factor  while  the  D value,  thermal 
diffusivity  of  food,  surface  heat  transfer  coefficient  and  container  geometry  had  minor 
influence  on  the  optimum  processing  temperatures. 


CHAPTER  1 
INTRODUCTION 

Thermal  processing  is  a popular  method  for  the  preservation  of  foods  It  produces 
versatile,  convenient  and  shelf  stable  products.  The  food  is  hermetically  sealed  in  a 
container  and  heated  at  the  proper  temperature  for  a sufficient  time  to  attain  commercial 
sterility.  Refrigeration  after  processing  is  usually  not  needed. 

The  healing  and  cooling  processes  are  critical  control  points  in  thermal  processing 
of  low-acid  foods.  All  commercial  processors  in  the  U.S.  and  those  in  other  countries  that 
export  such  foods  to  the  U.S.  are  required  to  tile  a process  schedule  for  all  thermally 
processed  products  to  ensure  safety.  Generally,  the  target  organism  is  the  Clostridium 
boltdiiium  which  produces  deadly  toxin. 

The  calculation  of  the  thermal  process  time  to  achieve  adequate  bacterial  lethality 
is  one  of  the  few  processes  in  the  food  industry  which  relies  on  a mathematical  model  to 
ensure  the  safety  of  the  products.  This  requires  knowledge  of  the  actual  transient 
temperature  history  of  the  food  during  the  process  by  either  predicting  it  with  a heat 
transfer  model  or  by  actual  measurement 

Heat  transfer  models  for  simulation  of  conduction  heated  foods  such  as  canned 
tishes  and  meats,  baby  foods,  pet  foods,  pumpkin,  and  squash  require  the  solution  of  a 
second  order  partial  differential  equation.  This  can  be  approached  by  two  methods: 


analytical  and  numerical  solutions  The  disadvantages  of  the  analytical  solutions  are  their 
limitations  in  requiring  the  material  to  conform  to  a regular  geometry,  its  initial 
temperature  must  be  uniform,  the  surrounding  temperature  must  be  constant,  and  thermal 
properties  arc  assumed  constant.  These  restrictions  can  be  overcome  by  the  use  of 
numerical  methods. 

The  most  common  numerical  method  used  in  the  analysis  of  heat  transfer  in  food 
processing  is  tile  finite  difference  approach.  Teixcira  ct  al  (1969)  first  introduced  the  use 
of  the  finite  difference  model  to  simulate  the  thermal  processing  of  conduction  healed 
foods  packed  in  a cylindrical  container.  As  computer  speed  and  capabilities  increased,  the 
model  has  gained  popularity  in  thermal  processing  design  and  analysis.  Today,  it  can  be 
used  to  predict  the  product  temperature  history  at  any  specified  location  within  the 
container  for  any  set  of  processing  conditions  and  container  size.  This  method  proved  to 
be  successful  in  predicting  product  temperature  both  in  constant  and  in  variable  retort 
temperature  conditions.  However,  the  finite  difference  model  is  only  applied  to  a few 
regular  geometries.  The  most  common  shape  is  the  finite  cylinder. 

The  recent  introduction  of  the  microwave  oven  for  home  food  preparation  has 
created  a consumer  demand  for  convenience  foods  in  microwaveablc  containers,  ready  to 
reheat  and  serve  in  the  package.  To  meet  this  need,  the  food  industry  has  introduced  shelf- 
stable  prepared  foods  in  new  containers  such  as  lunch  bowls  or  dinner  trays.  These 
containers  not  only  provide  the  convenience  for  serving  but  also  provide  an  alternate  way 
for  products  to  be  presented,  which  make  them  more  appealing  to  the  consumer. 
Developments  in  food  packaging  regarding  the  shape,  size  and  materials  of  packages  are 


very  dynamic.  Since  safety  of  the  thermally  processed  food  is  crucial,  thorough 
understanding  of  heat  transfer  behavior  in  any  new  container  is  needed.  This  prompted  a 
growing  interest  in  the  development  of  heat  transfer  models  for  these  new  shapes. 

The  container  shapes  of  interest  in  this  work  are  conical  and  finite  elliptical 
cylinder.  The  conical  container  is  commonly  used  to  pack  various  kinds  of  self-served 
foods  such  as  ice  cream,  yogurt,  cream  etc.  The  shape  could  be  varied  from  a pointed 
cone  to  an  almost  finite  cylinder.  The  heat  transfer  model  for  the  conical  shape  could  not 
be  found  in  the  literature. 

The  finite  elliptical  cylinder  is  commonly  used  to  pack  sardines  and  other  seafood. 
There  are  few  studies  of  heat  transfer  models  in  elliptically  shaped  containers  for  thermal 
processing. 

The  numerical  heat  transfer  model  could  also  be  used  to  predict  nutrient  or  quality 
retention  in  the  food  during  thermal  processing.  Conduction  heated  foods  change  their 
temperature  slowly  and  display  a non-uniform  temperature  distribution  during  healing  and 
cooling.  The  severity  of  healing  in  a container  varies  from  the  wall  to  its  geometrical 
center.  Because  of  this,  a high  temperature  short  time  (HTST)  process  will  not  result  in 
maximum  quality  retention.  A very  high  retort  temperature  will  cause  severe  thermal 
degradation  of  the  food  near  the  surface  before  the  center  of  the  container  has  reached  a 
desired  lethality.  On  the  other  extreme,  a low  retort  temperature  will  cause  a great  loss  of 
nutrient  since  long  processing  times  are  needed  to  obtain  commercial  sterility  Since  there 
are  infinite  time-temperature  combinations  which  give  the  same  bacterial  lethality,  the  use 
of  the  heat  transfer  model  together  with  calculations  for  nutrient  retention  predictions  will 


allow  the  processor  to  optimize  processing  conditions  for  quality  retention.  This  is  cost- 
and time-efficient  because  it  can  be  achieved  without  requiring  many  experimental  runs. 

The  climate  of  high  competition  in  food  marketing  and  the  increasing  demand  for 
high  quality  products  have  increased  the  demand  for  optimization  studies.  Processes 
producing  high  quality  products  will  attract  a larger  share  of  the  growing  health  and 
quality  conscious  consumer  market  A number  of  methods  have  been  developed  to 
maximize  quality  retention  in  thermal  processing  but  these  studies  are  limited  to  foods 
packed  in  finite  cylindrical,  rectangular,  and  spherical  containers. 

The  objectives  of  this  study  were  to  develop  heat  transfer  models  for  conduction 
heating  of  foods  packed  in  conically  and  elliptically  shaped  containers,  and  to  use  the 
developed  models  for  optimization  of  quality  retention  in  such  containers.  In  order  to 
fulfill  the  objectives,  the  research  was  conducted  in  3 phases: 

I Development  of  mathematical  models  of  conduction  heat  transfer  in  conically 
and  elliptically  shaped  container  in  response  to  constant  or  dynamic  surface  temperatures, 
for  predicting  the  content  temperature  at  any  given  location  in  the  container  and  at  any 
given  process  time. 

2.  Experimental  validation  of  the  models  by  comparing  model  temperature 
predictions  with  heat  penetration  data,  for  various  surface  temperature  conditions. 

3.  Use  of  the  developed  models  to  calculate  volume  average  nutrient  retention 
under  various  process  conditions,  with  the  goal  of  quality  optimization. 


CHAPTER  2 
LITERATURE  REVIEW 

Thermal  processing  of  foods  was  invented  by  Nicolas  Appert  in  1809.  The  method 
was  a combination  of  sealing  foods  in  a hermetic  container  and  heating  to  eliminate  or 
reduce  microorganisms  and  /or  enzymes  causing  deterioration  or  health  risk. 

Thermally  processed  foods  remain  popular.  They  provide  a safe,  convenient, 
nutritious  and  economical  alternative  to  fresh  and  frozen  foods.  Shelf  life  of  commercially 
sterilized  foods  is  at  least  two  years  at  room  temperature. 

Containers  for  pre-packaged  thermally  processed  foods  can  be  cylindrical, 
elliptical,  pear  or  rectangular  shaped  metal  containers,  glass  containers,  retortable 
pouches,  one-piece  plastic  trays,  bowls  or  cylindrical  plastic  cans,  etc. 

Before  heat  transfer  in  thermally  processed  foods  was  well  understood,  the  time- 
temperature  necessary  to  obtain  a safe  product  was  found  by  trial-and-error.  In  1920, 
Bigelow  and  Bali  developed  the  first  scientifically  based  method  for  the  calculation  of 
minimum  safe  sterilization  processes  for  sterilization  of  foods  in  cylindrical  shaped 
cansfLopez,  1987).  Ball  (1923)  proposed  a graphical  method  to  evaluate  the  degree  of 
sterility  of  food  based  on  the  temperature  history  at  the  center  of  the  can.  Thermal 
processing  techniques  have  since  been  well  developed,  making  thermal  processing  one  of 
the  safest  methods  for  preserving  foods.  Research  in  this  area  is  active  and  continuing. 


The  focus  now  are  maximizing  quality  and  retention  of  nutrient,  more  efficient  production, 
computerized  process  control,  more  convenience  and  appealing  packaging,  and  minimum 
energy  utilization. 

Heal  Transfer  Models  for  Conduction  Heated  Foods 

The  heat  transfer  model  for  conduction  heating  of  foods  requires  the  solution  of  a 
second-order  partial  differential  equation.  The  general  equation  for  unsteady-state  heat 
transfer  into  multi-dimensional  objects  is  described  as 

i(il)  . &!  . £l  , £1  (rectangular  geometry)  (2  1) 

a dt  r dr  dr 3 dy2  ^ 

Terms  are  explained  in  the  List  of  Symbols. 

The  solution  can  be  obtained  by  analytical  or  numerical  methods. 

Many  analytical  solutions  exist  for  conduction  heat  transfer  in  solids  of  different 
regular  geometries  such  as  cylinder,  sphere,  and  slab  (Carslaw  and  Jaeger,  1959).  Before 
the  computer  was  intensively  used,  a number  of  analytical  models  were  developed  to 
describe  the  temperature  history  at  the  slowest  heating  point  of  foods  in  regular  shaped 
containers  going  through  the  sterilization  process  (Gillespy,  1951;  Hicks,  1951). 


Hayakawa  (1969)  developed  a chart  using  an  analytical  method  to  facilitate  the  estimation 
of  the  healing  and  cooling  curves  at  the  cold  point  for  conduction  heated  foods  packed  in 
cylindrical  cans.  Experiments  were  conducted  to  validate  the  model  using  a metal  can 
heated  in  well  circulating  boiling  water  followed  by  cooling  the  can  in  well  circulating  cold 
water.  Infinite  heat  transfer  coefficient  conditions  were  assumed  for  both  heating  and 
cooling  stages.  Agreement  between  experimental  and  model  predicted  temperatures  was 
closer  during  the  heating  than  during  the  cooling.  This  was  explained  by  the  low  value  of 
the  heat  transfer  coefficient  during  cooling,  which  was  not  reflected  in  the  model. 

Analytical  solutions  have  limitations.  The  container  must  conform  to  a regular 
geometry,  initial  temperature  of  the  food  must  be  uniform,  surrounding  temperature  must 
be  constant,  heat  transfer  coefficient  and  the  thermal  properties  must  be  uniform  and 
constant  during  the  process. 

Mlffli.erifftl  Solutions 

Numerical  solutions  are  more  flexible  and  can  overcome  the  restrictive 
assumptions  of  the  analytical  solutions.  They  can  handle  problems  with  nonunifomi  initial 
temperature  distributions,  and  that  have  temperature-dependent  thermophysical  properties 
and  time-dependent  boundary  conditions.  The  numerical  simulation  of  thermal  processing 
of  foods  in  a cylindrical  can  was  first  introduced  by  Teixeira  ct  al.  ( 1 969).  They  used  a 
finite  difference  method  involving  Taylor's  expansions  to  find  the  solution  of  the  two- 
dimensional  partial  differential  equation  that  describes  conduction  heat  transfer  in  a 
cylindrical  can.  The  temperature  of  the  food  was  expressed  os  a function  of  the  outside 


temperature,  the  initial  product  temperature,  the  location  inside  the  container,  the  thermal 
diflusivity  of  the  product,  and  the  processing  time.  The  volume  was  subdivided  into  small 
volume  elements  (e.g.  concentric  rings  with  rectangular  cross  sections  for  a cylinder). 
Every  element  contained  a node  as  a reference  point.  By  assigning  boundary  and  initial 
conditions  to  every  node,  the  new  temperature  of  the  node  is  calculated  after  a short  time 
interval  (At),  and  the  procedure  is  repeated  until  the  end  of  process  time  (Teixeira,  1992). 

.Finnc-Ri.ncrcncc.  Method  with  Energy  Balance  Approach 

An  energy  balance  method  could  be  used  to  develop  finite  difference 
representations.  Instead  of  using  the  partial  differential  method  to  derive  finite  different 
equations,  an  energy  balance  equation  was  developed  for  every  volume  element.  The 
statement  of  the  energy  balance  with  no  energy  generation  is 

The  change  in  the  internal  The  net  heat  transfer  into  the 

energy  of  the  volume  element  = element  during  At 

during  the  time  interval  (At) 

The  typical  grid  pattern  for  a finite  cylinder  with  a convective  boundary  condition 
is  shown  in  Figure  2. 1 . The  temperature  at  each  node  represents  the  temperature  at  that 


location  in  the  volume  element. 


Figure  2.2  An  interior  node 


The  interior  node,  i,  j is  surrounded  by  4 adjacent  nodes(Figure  2.2).  The  overall 


The  energy  balance  equations  are  derived  from  the  physical  relationship  of  the 
elements.  It  is  clear  to  understand  and  provides  a meaningful  check  of  accuracy  for  the 
derived  differential  equations  and  the  resulting  computer  program  (Clausing.  1989).  The 
method  can  effectively  model  irregularly  shaped  objects,  temperature-dependent 
properties  and  phase  or  volume  change  in  heat  conduction  problems  (Chau  and  Snyder, 
1988;  Schreiber,  1990;  Sheen  and  Hayakawa,  1991;  Sheen  et.al.,1993;  Erdogdu,  1996  and 


Welt  et.al.,  1997). 


Heat  convection  occurs  whenever  a surface  t 


[ in  contact  with  fluid  at . 


different  temperature.  The  expression  used  to  determine  convection  heat  transfer  rate  is 


In  the  energy  balance  approach,  when  the  volume  elements  facing  the  convective 
boundary,  the  nodes  could  be  set  as  capacitance  surface  nodes(CSN)  or  non-capacitance 
surface  nodes  (NCSN)  (Clausing,  1981).  The  capacitance  surface  nodes  are  nodes  on  the 
surface  of  volume  elements  with  thermal  mass.  The  temperatures  of  the  elements  are 
changed  instantly  with  changing  outside  temperature,  which  could  result  in  error  especially 
at  small  values  of  time.  A small  element  is  required  to  reduce  such  error,  resulting  in 
decreasing  maximum  time  step  and  increasing  computational  time. 


Figure  2.3  Node(ij)  with  non  capacitance  surface  nodes  on  its  right  and  bottom  surfaces 


C ■ * a ( r-rj 


(2.4| 


The  non  capacitance  surface  nodes  are  the  nodes  on  the  surface  of  the  volume 
element  that  had  no  thermal  mass.  In  Figure  2.3,  node  (ij)  has  NCSN  on  the  right  and  the 


bottom  surface.  Heat  transfers  through  the  surface  to  the  node  (i  j)  is  analyzed  as  a 
combination  of  resistance  to  heat  transfer  involving  both  convection  and  conduction  as 


By  setting  the  NCSN,  accurate  results  are  obtained  and  the  time  step  restriction 
under  high  convective  heat  transfer  coefficient  condition  is  eliminated.  The  works  involved 
with  NCSN  include  Bellagha  and  Chau(l985),  Chau  and  Gaffiicy(l990),  Silva  et  al  (1992 
a).  Erdogdu  (1996),  and  Welt  ct  al.(!997). 

Convective  Heat  Transfer  Coefficient 

The  convection  heat  transfer  coefficient  (h)  is  a complex  function  of  geometry, 
fluid  flow  and  fluid  properties.  An  infinite  convective  heat  transfer  coefficient  with 
negligible  surface  resistance  is  a valid  assumption  when  condensing  steam  is  used  to  heat 
foods  in  metal  containers.  By  this  assumption,  temperature  of  the  food  at  the  surface  is 
equal  to  the  temperature  of  the  heating  medium.  When  prepacked  foods  are  cooled  with 
water  or  heated  with  a steam/air  mixture,  a finite  surface  heat  transfer  coefficient  is 
assumed  (Tucker  and  Clark,  1990;  Bhowmik  and  Tandon,  1987;  Lebowilz  and  Bhowmik, 
1989,1990;  Silva  etal.,  1992b). 

When  foods  are  packed  in  a high  heal  barrier  material,  e g.,  glass,  or  plastic 
containers,  the  overall  heat  transfer  coefficient  (U)  is  used.  This  value  is  the  combination 


of  outside  surface  heat  transfer  coefficient  and  the  thermal  resistance  of  the  packaging 
material  (Lebowitz  and  Bhowmik,  1990).  Sheen  et  al.  (1993)  used  a U value  of 
700  W / m 1 K in  the  finite  volume  model  for  simulating  heat  transfer  in  apple  sauce  packed 
in  a bowl  shaped  plastic  container  (5.3  mm  thickness)  processed  in  an  autoclave,  Lebowitz 
and  Bhowmik  (1990)  determined  overall  heat  transfer  coefficients  for  the  come-up, 
heating  and  cooling  cycles  of  thermal  processing  for  retortablc  pouches  heated  by 
circulating  hot  water  under  overriding  air  pressure  They  found  that  U values  for  these 
different  processing  cycles  were  not  significantly  different  and  were  within  the  range  of 
174-243  W / m!  K.  Since  the  pouch  materials  are  very  thin,  they  also  found  that  U values 
for  either  aluminum/fbil/polymer  laminate  or  all  plastic  pouches  were  not  significantly 
different. 


There  are  some  research  on  the  application  of  finite  difference  methods  to  solve 
transient  heat  transfer  problems  involving  irregularly  shaped  geometries.  Manson  et  al. 


(1974)  developed  a 3-dimensional  finite  difference  model  for  conduction  heated  foods  in  a 
pear  shaped  container.  However  any  attempt  to  validate  the  model  was  not  mentioned. 
Simpson  et  al.  (1989)  developed  a 3-dimensional  finite  difference  model  to  evaluate 
thermal  processes  for  conduction  heated  foods  in  an  oval  shaped  container.  They  used 
Equation  2. 1 in  rectangular  coordinates,  and  the  volume  was  subdivided  into  small 
rectangular  elements.  To  make  the  grid  fit  the  curved  boundaries,  a proportionality  factor 
was  applied  to  the  nodes  near  the  surface  of  the  oval.  In  this  study,  the  predicted 


temperature  within  the  oval  shape  was  validated  with  an  analytical  solution,  but  not 
experimentally.  Chau  and  Snyder  (1 9S8)  developed  a Unite  different  model  for  thermally 
processed  shrimp  using  the  energy  balance  technique  with  non  capacitance  surface  node 
A shrimp  was  approximated  by  a short  cylinder  followed  by  a right  circular  cone.  A stair 
step  approach  was  used  to  approximate  the  radial  boundary  of  the  cone.  The  volume 
elements  were  series  of  concentric  rings  stacked  upon  each  other  with  each  consecutive 
set  of  rings  containing  one  less  node.  The  experimental  validations  of  a center  temperature 
of  a shrimp  during  heating  and  cooling  were  conducted  in  a water  bath.  A good  agreement 
of  experimental  and  simulated  results  was  found.  Sheen  et  al.  (1993)  applied  finite  volume 
method  to  a curvilinear  boundary  to  develop  a model  for  conduction  heated  food  in  a 
bowl-shaped  plastic  container  An  experimental  validation  was  carried  out  in  an  autoclave 
using  bentonite  and  apple  sauce.  They  found  satisfactory  agreement  between  the  model 
prediction  and  experimentally  measured  temperatures  at  two  points  within  the  bowl  of 
bentonite  solution.  Erdogdu(  1 996)  developed  a finite  difference  model  using  the  energy 
balance  technique  for  heat  transfer  in  an  elliptical  cylinder  with  the  method  of  heat 
flowlines  as  suggested  by  Eshleman(1976)  and  Pham(l99l).  There  was  close  agreement 
between  experimental  and  simulated  temperature  histories  at  the  center  and  an  off-center 
point  of  an  elliptical  cylinder.  He  found  that  it  is  also  possible  to  accurately  predict  center 
temperatures  of  an  elliptical  cylinder  with  a right  circular  cylinder  of  an  equal  surface-area- 
to-volume  ratio  when  the  eccentricity  of  the  cross  section  ellipse  < 0.7. 

With  increasing  computer  speed  and  capability,  numerical  methods  have  gained 
more  popularity  and  have  been  widely  used  by  several  researchers  to  predict  temperature 


distribution  of  conduction  heated  foods.  Table  2. 1 shows  some  uses  of  the  Unite 
difference  method  in  thermal  processing  of  foods. 

A finite  element  method  could  also  predict  temperature  distribution  of  foods 
during  thermal  processing.  Although  the  finite  element  method  is  more  flexible  for 
irregular  containers  than  finite  difference,  it  requires  more  computer  time,  and  is  more 
difficult  to  set  up.  Tendon  and  Bhowmik  (1986)  developed  a finite  difference  model  for 
predicting  temperature  for  conduction  heated  foods  packaged  in  retortable  pouches  where 
the  actual  shape  was  considered.  A proportionality  factor  was  applied  to  the  elements 
near  the  edge  of  the  container.  They  also  used  the  finite  element  method  to  model  heal 
transfer  in  the  same  object.  Temperature  prediction  validation  was  conducted 
experimentally  with  a bentonite  solution  in  a retort.  They  found  that  with  the  same 
accuracy,  the  finite  difference  method  required  4 times  less  computer  time  and  was  easier 
to  set  up. 

Application  of  Numerical  Heat  Transfer  Model  in  Thermal  Processing  of  Foods 

Several  applications  of  the  numerical  model  to  simulate  temperature  history  at  a 
given  location  of  packed  foods  undergoing  thermal  processing  are  discussed  by  Datta  et 
al  ( 1986)  and  Teixeira  (1992).  A model  could  predict  product  temperature  either  when  the 
retort  temperature  was  constant  or  varied  in  any  prescribed  manner  throughout  the 
process.  The  use  of  these  models  has  become  invaluable  since  it  is  not  necessary  to  carry 
out  repeated  heat  penetration  tests.  These  models  could  be  used  for  simulating  the  process 
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conditions  experienced  in  continuous  sterilizer  systems,  in  which  containers  pass  front  one 
chamber  to  another,  experiencing  a changing  retort  temperature  at  the  can  wall.  The 
models  could  also  be  used  to  rapidly  evaluate  an  unscheduled  process  deviation,  such  as 
when  an  unexpected  drop  in  retort  temperature  occurs  during  processing,  by  predicting 
product  center  temperature  in  response  to  that  variation,  and  evaluate  a sterilizing  value 
based  on  that  prediction. 

Another  essential  use  of  the  models  is  to  find  the  optimal  process  condition  giving 
the  maximum  nutrient  or  quality  retention  while  maintaining  a desired  safety  of  the 
product.  The  model  for  this  application  should  be  able  to  predict  the  temperature  for  all  of 
the  locations  within  the  volume  throughout  the  process,  and  use  these  data  to  calculate 
the  retention  of  quality  or  nutrient.  Without  the  use  of  the  models,  optimization  of  nutrient 
retention  becomes  almost  impossible  since  it  would  require  a series  of  heat  penetration 
tests  and  experiments  to  measure  the  retention  of  the  quality  attributes  of  interest  after  the 

Effect  of  Heat  pn  Thermo-labile  Substances  in  Food 

Thermo-labile  factors  of  foods  that  degrade  according  to  a first  order  kinetics 
include  most  microorganisms  and  their  spores,  enzymes,  vitamins  and  other  sensory 
quality  factors  (e.g.  color,  texture  and  flavor).  The  general  equation  to  describe  the  first 
order  degradation  kinetics  can  be  written  as 


In  order  to  define  the  temperature  sensitivity  of  the  kinetic  parameters,  two 
methods  are  used:  Arrhenius  and  thermal  death  time  (TDT)  methods.  The  Arrhenius 
method  is  based  on  a thermodynamic  approach.  The  natural  logarithm  of  reaction  rate 
constant  k is  related  to  the  reciprocal  of  the  absolute  temperature,  with  activation  energy 
(Ea)  found  from  the  slope  of  the  fitted  line.  The  method  has  been  widely  used  to  describe 
chemical  kinetics. 

The  TDT  concept,  first  introduced  by  Bigelow  et  al.(I920),  is  based  on  an 
empirical  approach  to  the  temperature  dependence  of  the  first  order  destruction  rate  of 
microorganisms.  This  method  has  been  subsequently  used  in  numerous  studies  leading  to 
its  widespread  application  in  thermal  process  calculations.  In  the  TDT  method,  the  D 
value  is  defined  as  the  time,  at  a constant  temperature,  necessary  to  reduce  90%  of  the 
initial  concentration.  The  relationship  between  D and  k is 


(2.8| 


Temperature  of  121 . 1°  C is  commonly  used  os  a reference  temperature  in  thermal 
processing.  The  D value  at  this  temperature  is  called  D,,,  A temperature  change  necessary 
to  cause  a 90%  change  in  D value  is  defined  as  the  Z value.  Reactions  having  a small  Z 


value  are  highly  temperature  dependent:  changing  temperature  has  a great  effect  on 


changing  the  reaction  rate.  On  the  other  hand,  reactions  with  high  Z values  have  low 
thermal  dependency  and  reaction  rates  are  less  sensitive  to  temperature  changes. 

The  values  for  D,[(  and  Z of  some  microorganisms,  spores  and  quality  attributes 
are  shown  in  Table  2.2.  The  values  for  a given  nutrient  or  quality  factor  depend  on 
intrinsic  properties  of  the  food  (e.g.  composition,  pH,  water  activity,  oxidation-reduction 
potential),  extrinsic  factors  (e.g.  gas  composition  inside  the  container)  and  their 
interactions  (Silva  et  a!.,  1992a).  1)1656  parameters  must  be  used  in  the  same  environment 
and  temperature  range  in  which  the  kinetic  data  were  obtained. 

E I.  on  of  Safety  in  Thermally  Processed  Foods 

The  main  objective  of  thermal  processing  is  to  destroy  spoilage  and  pathogenic 
microorganisms  together  with  their  spores  or  toxins  in  order  to  obtain  safe  and  shelf-stable 
products.  The  amount  of  heating  required  depends  on  the  resistance  to  heat  of  the 
microorganisms  which  may  be  present  in  the  container,  on  the  rate  of  heat  penetration  into 
the  "cold  point"  of  the  container,  and  on  the  desired  level  of  lethality  required. 

Concept  of  Cold  Point  Lethality 

Bigelow  et  al.(l920)  first  proposed  the  general  method  to  evaluate  the  lethality  of 
bacterial  spores  during  thermal  processing.  The  calculation  of  lethality  is  based  on  actual 
product  temperature  at  the  slowest  heating  point  (cold  point)  and  integrates  the  total  lethal 
effect  of  time-temperature  combinations  for  the  entire  process. 

Total  lethality  is  described  mathematically  by  the  F0  value.  It  is  defined  as  time  at 


the  reference  temperature  (generally  1 2 1 . 1 °C ), ; 


i attainment  of  this 


temperature  and  instantaneous  cooling  afterwards,  that  would  give  the  same  degree  of 
sterility  as  the  process.  F„  value  can  be  expressed  as 


Clostridium  botulinum  is  the  most  important  target  spore  forming  microorganism 
in  relation  to  the  safety  of  low  acid  heat-processed  foods.  A process  would  be  described  as 
safe  when  it  reduces  the  probability  of  the  survival  of  this  organism  to  a very  small  value, 
such  as  not  more  than  one  in  a trillion  cans  (I O'1)  (Gillespy.  1951).  Pflug  and  Odlaug 
(1978)  reviewed  the  literature  and  concluded  that  Z - I0°C,  D„,  = 0.2  min  could  be  used 
in  thermal  process  calculations  to  achieve  F„  = 3 min.  This  process  ensures  the  safety  of 
low-acid  canned  foods  from  botulism.  F„  for  some  current  commercial  processes  are 
presented  in  Table  2.3. 

Ball  (1923)  derived  a formula  for  a heat  penetration  curve  at  the  center  of  a 
cylinder  based  on  the  analytical  solution  of  conduction  heat  transfer.  This  was  called  the 
“formula  method."  Even  though  the  method  is  easier  to  perform  manually,  it  has  many 
restrictive  assumptions,  e.g.,  the  heal  penetration  curve  has  to  be  a straight  line,  the 
heating  and  cooling  temperature  must  remain  constant. 
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> : Lopez  (1987) 


With  the  progresses  that  have  evolved  in  computer  technology,  the  evaluation  of 
target  F„  by  the  general  method  has  become  more  popular  due  to  its  versatility.  Computer 
programs  available  commercially  require  only  time-temperature  data  at  the  cold  point  to 
evaluate  target  lethality  automatically. 

Position  pf  the  Cold  Point 

The  position  of  the  cold  point  of  food  undergoing  thermal  processing  is  crucial  tor 
lethality  calculation.  The  F„  value  at  the  cold  point  is  attributed  to  the  whole  container,  on 
the  principle  that  a process  which  is  safe  for  the  point  of  slowest  heating  must  be  safe  for 
the  whole  container.  Inaccuracy  in  locating  the  position  of  the  cold  point  may  result  in 
under-processing.  The  position  of  cold  point  for  convection  and  conduction  heated  foods 
is  different.  Agitation  or  rotation  of  a container  during  processing  also  affects  the  cold 
point  position. 

The  cold  point  for  conduction  heated  foods  with  no  agitation,  packed  in 
symmetrica]  shaped  containers  made  from  one  uniform  material,  e g.,  cylinder,  rectangular 
or  elliptically  shaped  metal  can  as  well  as  aluminum  retortable  pouch,  is  at  the  geometrical 
center.  This  is  not  true  for  asymmetrical  shapes  such  as  pear,  bowl  or  cone,  or  a container 
made  of  more  than  one  material,  such  as,  a plastic  can  with  a metal  lid.  Sheen  et  al. 

(1993)  found  that  the  cold  point  of  a bowl-shaped  container  with  a metal  lid  (top  radius : 
42  mm,  bottom  radius  = 28  mm,  height  = 57  mm,  with  curved  boundary)  during  heating  in 
a circulated  superheated  water-type  retort  was  located  at  a fraction  of  0.46  of  the  height 
of  the  container  from  the  bottom,  on  the  rotation  axis.  The  cold  point  moved  slightly 


toward  the  bottom  during  cooling  with  water  due  to  changing  heat  transfer  coefficient. 
Berry  and  Bush  (1988)  reported  that  the  cold  point  of  conduction  heated  food  in  a 
cylindrically  shaped  plastic  container  with  the  metal  lid  up  during  processing  was  slightly 
below  the  center  of  the  can. 

Processing  conditions  and  thcrmophysical  properties  of  foods  also  affect  the 
position  of  the  cold  point.  Naveh  et  al.  ( 1 983)  developed  a conduction  heat  transfer  model 
for  a product  packed  in  a glass  jar  with  head  space.  They  found  the  products  of  high 
thermal  diffiisivity  had  a slight  shift  in  the  location  of  the  cold  point  towards  the  head 

Concept  of  Volume  Average,  lethality 

During  a thermal  process,  the  severity  of  a thermal  process  in  a conduction  healed 
food  decreases  going  from  the  geometrical  center  to  the  wall  of  the  container.  Slumbo 
(1953)  proposed  a method  to  evaluate  volume-average  bacterial  spore  retention. 

He  described  this  gradient  of  severity  by  picturing  a series  of  imaginary  nested 
containers  inside  the  real  container.  These  containers  were  considered  to  be  cylindrical  in 
shape  and  decreasing  in  size  from  the  wall  to  the  geometrical  center  of  the  real  container, 
and  the  severity  of  heating  was  the  same  at  all  points  in  these  containers.  He  called  them 
iso-j  regions.  The  temperature  history  during  thermal  processing  for  one  point  of 
each  iso-j  region  was  calculated  by  an  analytical  method.  The  bacterial  spore 
retention  corresponding  to  that  point  or  its  iso-j-  region  was  evaluated  by  the 
general  method  for  thermal  process  calculation.  The  spore  retention  was  then 


containers  is  very  difficult(Pfiug  and  Odlaug.  1978).  However,  the  concept  became 
valuable  since  it  could  be  applied  to  determine  the  retention  of  any  heat  labile  component 
in  food,  e.g.  nutrients,  sensory  qualities,  which  follow  first  order  rate  of  destruction. 


Volume  Averarte  Quality  Retention 

Ball  and  Olson  (I9S7),  Jen  et  al  (1971),  and  Barreiro  et  al. (1984a)  extended  the 
volume  average  concept  to  develop  manual  methods  for  estimating  nutrient  degradation  in 
conduction  heated  foods  in  cylindrical  and  rectangular  containers.  The  methods  require 
constant  processing  temperature,  and  the  linear  heat  penetration  curve  is  assumed. 

When  Teixcra  et  al  (1909)  introduced  a finite  difierencc  model  to  predict 
temperature  at  any  location  al  any  given  time  in  the  containers  during  thermal  processing, 
they  combined  the  model  with  the  volume  average  concept  to  estimate  bacterial  spore  and 
thiamine  retention  afier  processing.  By  dividing  the  overall  volume  into  small  volume 
elements,  the  retention  afier  a small  time  interval  At  can  be  written  as 


where  N,  and  N ,.A,  are  the  number  or  concentration  at  time  t.  and  after  a short 
time  interval  (At),  respectively.  Dr  is  a decimal  reduction  value  which  is  temperature- 
dependent  and  will  be  updated  with  temperature  change  as 
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Optimization  in  Thermal  Processing 


The  optimization  in  thermal  processing  requires  the  identification  of  a process 
variable  that  can  be  controlled,  the  requirements  that  must  be  met  (the  constraint),  and  an 
objective  function  (Norback,  1980).  The  constraint  for  optimizing  the  heat  sterilization 
process  is  the  target  lethality  value  (F0)  which  is  required  to  obtain  a safe  product 

Several  objective  functions  have  been  considered  in  attempts  to  optimize  heat 
sterilization  processes.  Barreiro  et  al. (1984b)  determined  an  optimal  temperature  to 
process  canned  pea  puree  while  minimizing  energy  utilization  in  a batch  retort,  while 
maintaining  the  same  lethality  value.  However,  foods  that  are  processed  based  on 
minimum  energy  utilization  may  result  in  poor,  unacceptable  quality.  Therefore,  the  aim  of 
thermal  process  optimization  is  moving  toward  obtaining  maximum  nutrient  and/or  quality 
retention  (Hendrickx  et  al.,  1993). 

Microorganisms  and  their  spores  arc  much  more  heat  sensitive  than  nutrients  and 
sensory  factors  (Lund,  1977).  In  liquid  or  semi-liquid  products  which  have  uniform 
temperature  distribution  during  healing  and  cooling,  the  High  Temperature  Short  Time 
(HTST)  process  or  even  the  more  extreme  Ultra  High  Temperature  (UHT)  process  results 
in  higher  nutrient  retention.  However,  these  processes  could  not  be  applied  to  conduction 
heated  foods  (e.g.  tuna,  ham,  salmon,  pork  and  beans,  corned  beef  and  canned  pumpkin) 
due  to  non-uniform  temperature  distribution.  Very  high  temperatures  will  cause  severe 
thermal  degradation  of  quality/nutrients  near  the  container  wall  long  before  the  food  at  the 
center  of  the  container  has  risen  in  temperature.  On  the  other  extreme,  a relatively  low 


[ temperature  will  also  cause  great  quality  loss  because  of  the  long  times  necessary  to 


achieve  sterility  (Thijssen  et  al„  1978), 

There  is  an  infinite  number  of  time-  temperature  combinations  to  achieve  a desired 
lethality.  These  combinations  result  in  different  nutrient  and  quality  retention.  Therefore, 
optimization  of  processing  conditions  is  needed  to  obtain  good  quality  products. 

Choice  of  an  Objective  Function  to  Maximize  Quality  Retention 

Two  objective  functions  are  commonly  used  to  obtain  optimal  thermal  process 
conditions  to  maximize  quality  retention:  maximizing  volume  average  retention,  or 
minimizing  cook  value.  These  could  be  used  for  overall  volume  ( N/Numl  or  and/or 
surface  ( A'A',,,,,  or  CIUJ  of  the  container.  Table  2.4  presents  different  objective  functions 


Container  Shape,  Size  and  Material 

Teixeira  et  al.  (1975)  used  a computer  model  to  predict  thiamine  retention  in 
different  height-to-diameter  ratios  and  found  that  retentions  greater  than  60%  were 
achieved  for  cans  with  height-to-diamctcr  ratio  less  than  0.2  and  those  with  ratios  greater 
than  10.0.  A minimum  retention  of  about  40%  oocurred  with  ratios  around  I (standard 
size  cylindrical  cans). 

In  flat  containers  such  as  rctonable  pouches  or  trays  heat  will  transfer  faster  than 
in  a cylindrical  can  with  equivalent  volume.  Higher  optimum  sterilization  temperatures  and 
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belief  surface  and  average  quality  retention  arc  achieved  (Ohlsson,  1980a  ; Silva  et  al. 
1994a). 

Packaging  materials  have  an  important  role  in  the  rate  of  heat  transfer  into  the 
containers  and  therefore  nutrient  retention  after  processing.  Berry  and  Bush  (1988)  and  Lu 
et  al.  (1991)  studied  the  rate  of  heat  penetration  into  two  piece,  cylindrical^'  shaped  plastic 
containers  with  metal  lids  in  comparison  to  metal  cans  of  equal  dimensions.  They  found 
that  the  rate  of  heat  penetration  into  the  plastic  containers  was  lower  and  quality  retention 
was  poorer. 

The  thermal  resistance  or  Z value  has  a significant  effect  on  optimal  processing 
temperature  when  trying  to  maximize  average  retention  or  minimize  average  cook  value  of 
quality  factors  in  conduction  heated  foods.  Quality  factors  with  high  Z values  require  higlt- 
tcmpcraturc-short-time  processes.  This  is  in  contrast  to  those  of  lower  Z values  where  low 
temperature  long  time  process  is  preferred  (Teixeira  et  al.,  1969;  Ohlsson,  1980b; 

Silva  et  al.,  1992a) 

D-  value 

Teixeira  et  al.(l  969)  studied  the  effect  of  varying  the  D„,  value  from  146.6  to 
188.37  to  202.3  min  on  the  optimal  temperature  of  pea  puree  and  concluded  that  this 
variable  was  not  significant.  They  found  that  the  higher  D„r  value  simply  increased  the 
percent  retention.  However,  they  coasidered  only  a narrow  range  of  Dre(,  and  this  range 
does  not  include  the  D,,,  values  for  sensory  quality  factors  (10  to  SO  min).  Silva  et  al 


( 1 992a)  extended  die  range  of  Dn(  values  (between  5 and  1 000  minutes)  and  they  found 
that  D",  less  than  I SO  min  has  a significant  effect  on  optimal  temperature.  The  quality 
factors  with  low  DM  prefer  high  optimal  temperature-short  time  processes  They  also 
stated  that  the  use  of  volume  average  cook-value  which  does  not  take  into  account  D„, 
may  result  in  underestimation  of  optimal  processing  temperature. 


Most  optimal  sterilization  processes  have  been  determined  manually  (Teixeira  et 
al„  1969,  1975;  Ohlsson,  1980a;  Thijssen  and  Kochen,  l980;TukerandHoldsworth, 
1990;  Hendrickx  et  al.,  1992).  The  final  quality  retention  is  calculated  for  different 
constant  temperatures  or  tompcraiurc  profiles,  and  plotted  as  a function  of  the  design 
variable.  The  optimum  condition  is  obtained  by  graphical  observation  (Silva  et  al.,  1993). 

In  1979,  Saguy  and  Karel  were  the  first  to  introduce  a formal  optimization  theory 
to  obtain  maximum  quality  retention.  They  applied  the  continuous  maximum  principle  of 
Pontryagin  and  a multi-iterative  method  to  determine  the  best  time-temperature  profile  in 
order  to  optimize  thiamine  retention  in  thermal  processing.  f\adkarni  and  Hatton  ( 1 98S) 
used  the  minimum  principle  of  optimal  control  theory  with  the  goal  of  optimizing  the 
nutrient  retention  for  a given  reduction  of  microorganisms.  Silva  et  al.  (1992b)  used  the 
Davis,  Swann  and  Campey  method  to  calculate  optimal  constant  temperatures. 

Hendrickx  et  al.  (1990)  used  an  empirical  approach  to  calculate  optimal 
temperatures  for  maximizing  surface  cook  value  with  infinite  surface  heat  transfer 
coefficients  for  one  dimensional  shapes  (infinite  slab,  infinite  cylinder  and  sphere).  Using 


traditional  regression  analysis,  they  developed  empirical  equations  to  relate  optimal 
temperatures  to  various  product  properties  (thermal  diflusivity  and  Z-value  of  quality 
factors),  processing  conditions  (geometry  and  dimensions  of  the  food,  surface  heat 
transfer  coefficient,  initial  product  temperature  and  retort  come  up  time)  and  processing 
criteria  (target  F„  value).  Hendrickx  et  al.  ( 1 993)  extended  the  correlations  for  more 
generalized  conditions  accounting  for  cooling  lethality  as  well  as  retort  come  up  lime. 
They  found  that  heating  medium  come  up  time  and  the  initial  product  temperature  have  a 
minor  effect  on  the  optimal  processing  temperature.  Silva  et  al.(!992b)  presented  the 
correlations  for  optimal  sterilization  temperature  for  conduction  heated  foods  with  finite 
surface  heat  transfer  coefficients.  Sablani  (1995)  applied  a neural  network  approach  to 
compute  optimal  sterilization  temperature  and  their  corresponding  processing  time  for 
finite  cylinder  cans. 

Experimental  Validation  of  Quality  Retention 

Some  experiments  have  been  conducted  to  validate  the  quality  prediction  models 
The  validation  requires  the  selection  of  a quality  parameter,  the  kinetic  data  for  this 
parameter,  heat  penetration  characteristics  in  a filled  container  and  processing  using 
different  sterilization  temperatures  ( Silva  et  al.,l994b). 

The  selected  quality  parameters  for  validation  tests  included  surface  lightness  of 
tuna  (Banga.1993),  Hunter  color  scale  values  of  pea  puree  (Shin  and  Bhowmik,  1995), 
and  acid  catalyzed  sucrose  hydrolysis  in  bentonite  solution  (Silva  et  al.,l994b).  The 
accuracy  of 
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ihe  predictions  relies  on  the  accuracy  of  predicting  the  temperature  at  each  and  every 
position  during  the  heating  and  cooling  process,  and  on  the  accuracy  of  the  kinetic  data. 

The  objectives  of  this  study  were  to  develop  conduction  heat  transfer  models  for 
elliptical  and  conical  shaped  containers,  and  apply  the  models  for  thermal  processing  to 
estimate  quality  retention.  To  achieve  the  objectives,  the  steps  taken  were  to  develop  the 
models  using  the  finite  difference  method  based  on  the  energy  balance  approach;  to 
evaluate  the  accuracy  of  the  models,  the  temperature  distributions  within  a cylinder  were 
simulated  with  the  developed  finite  elliptical  and  conical  models  and  the  results  compared 
to  those  obtained  from  the  validated  numerical  finite  cylindrical  model.  For  the 
experimental  validation,  acrylic  rods  cut  into  various  shapes  and  sizes  of  cone  frustum  and 
finite  elliptical  cylinder  were  fitted  with  thermocouples  to  measure  interior  temperatures  at 
various  locations.  The  experiments  were  carried  out  under  constant  and  variable  outside 
temperature  environments.  Experimental  temperatures  at  different  locations  were 
compared  to  simulated  temperatures. 

The  validated  models  were  applied  to  calculate  the  volume  integrated  quality 
retention.  The  influences  of  the  surface  heat  transfer  coefficient,  the  thermal  difiusivity  of 
food,  and  the  different  container  geometry  and  size  on  the  retention  of  quality  factors  with 
different  kinetic  parameters!  D„,  tmd  Z)  values  were  investigated.  The  goal  of  the 
investigation  was  to  optimize  the  quality  of  the  food. 


CHAPTER  3 

MATERIALS  AND  METHODS 
Development  of  ihc  Numerical  Heal  Transfer  Models 

The  finite  difference  method,  based  on  energy  balance  approach,  was  used  to 
develop  conduction  heat  transfer  models  in  foods  contained  in  conically  and  clhpucally 
shaped  containers. 

Basic  assumptions  of  the  models  include 

• Heat  transfer  into  the  food  is  by  pure  conduction 

• The  food  product  is  homogeneous  and  isotropic,  i.e.,  its  properties  are 
constant  with  location  and  temperature. 

• The  initial  temperature  of  the  food  is  uniform. 

• The  surface  heat  transfer  coefficient  can  be  finite  or  infinite. 

The  first  step  in  developing  the  heat  transfer  model  is  to  divide  the  overall  volume 
into  small  volume  elements.  A node  is  assigned  to  every  volume  element  as  a reference 
point  and  is  normally  located  at  the  center  of  the  element.  The  temperature  at  the  node 
location  is  the  temperature  at  that  location  in  the  volume  element.  Each  element  has 
imaginary  boundaries  around  the  node,  and  the  resulting  grid  forms  a continuous  network 
throughout  the  volume.  The  geometrical  properties  of  the  elements  for  the  entire  region 


are  required  to  develop  energy  balance  equations.  These  are:  distance  between  adjacent 
nodes,  the  volume  of  the  element  surrounding  each  node,  and  the  area  of  heat  transfer  for 
the  node.  The  following  properties  arc  also  needed:  thermal  conductivity,  density,  and 
heat  capacity  of  the  food,  surface  heat  transfer  coefficient,  and  outside  temperature  The 
next  step  is  to  develop  an  energy  balance  equation  for  every  element.  It  was  assumed  that 
heat  was  transferred  perpendicularly  into  a node  from  2 directions:  radial  and  axial. 
According  to  the  definition  of  Clausing  ( 1 987),  the  interior  nodes  are  defined  as 
capacitance  nodes.  They  have  a volume  and  heat  capacity.  The  nodes  at  the  surface  may 
be  capacitance  or  noil-capacitance  nodes.  The  nodes  with  similar  geometry  and  boundary 
condition  are  grouped  together  to  facilitate  computer  coding. 

In  this  study,  an  explicit  solution  was  used  to  calculate  the  nodal  temperatures 
using  the  energy  balance  equations.  The  temperature  after  a small  time  step  depends  on 
the  temperatures  before  this  time  step.  The  maximum  time  step  is  limited  by  the  second 
law  of  thermodynamics.  Therefore,  it  is  necessary  to  develop  the  equations  for 
calculating  the  maximum  time  step. 

Qrid  generation  in  a Cong  Frustum 

The  dimensions  of  a cone  frustum  arc  given  as:  big  radius  (R),  small  radius  (r)  and 
height  (H).  The  cone  is  subdivided  into  a series  of  concentric  rings  stacked  on  top  of  each 
other(Figure  3 .1).  The  volume  elements  at  the  side  surfaces  have  a shape  of  triangular 
cross-section  ring.  The  number  of  volume  elements  increase  by  one  for  every  axial  level. 
Because  of  the  axial  symmetry,  only  one-half  of  it  will  be  considered.  As  shown  in  Figure 
3.2,  nodes  are  at  the  center  of  every  volume  element  as  a reference  point,  except  for  the 
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Figure  3. 1 Grid  generalic 


: frustum 


Figure  3.2  Volume  elcmenl  definitions  in  a cone  frustum 


# Capacitance  nodes 
O Non  capacitance  surface  nodes 


nodes  at  the  surface.  The  interior  nodes  are  capacitance  nodes.  The  nodes  on  the  top  and 
bottom  surfaces  are  non-capacitance  surface  nodes,  while  nodes  on  the  side  surfaces  of  the 
cone  are  capacitance  surface  nodes.  The  number  of  nodes  in  the  R-  r region  equals  to 
number  of  nodes  in  height  (not  including  non-capacitance  surface  nodes).  The  expressions 
of  distances  between  nodes,  areas  of  heat  transfer  and  nodal  volume  are  presented  in 
Appendix  A. 


Ten  different  types  of  nodes  were  generated  within  a cone  frustum  as  shown  in 
Figure  3.2  Energy  balance  equations  and  the  calculated  maximum  time  steps  for  the 
volume  elements  were  as  follows: 

Node  Type  I 


Specification 

This  node  is  a non-capacitance  surface  node  (NCSN),  having  no  volume.  It  is 
located  at  the  top  surfacefFigure  3,3).  Heat  convects  from  the  top  surface  boundary  into 
the  node  and  conducts  to  the  node  at  the  bottom  Since  heat  can  not  be  stored  in  the 
node,  its  energy  balance  equation  is  written  as 


Figure  3.3  Representation  of  node  type 


(3.1) 


The  temperature  of  the  node  is 


(3.2) 


There  is  no  need  to  calculate  a maximum  time  step  for  this  type  of  node. 


This  node  is  a NCSN  similar  to  node  type  I except  it  is  located  at  the  bottom 
surface  (Figure  3.4).  The  temperature  of  the  node  can  be  written  as 


(3.3) 


Node  Type  3 


Figure  3.5  Representation  of  node  type  3 


This  node  is  a capacitance  node  located  at  the  top  surface.  The  top  of  the  node  is  a 
NCSN  with  convection  boundary.  It  is  surrounded  by  three  capacitance  nodes 
(Figure  3.5).  The  energy  balance  equation  for  this  node  is 


(3.4) 


It  is  possible  to  substitute  T' .....  in  equation  3.4  with  equation  3.2.  Since  the 
reference  node  is  now  node  type  3,  the  equation  3.2  is  then  changed  to 


(3.5) 


After  substituting  thermal  difTusivitv  ( a = ) and  solving  for  the  nodal  temperature 

P Cp 

in  the  next  time  step  ( T-''1 ) the  corresponding  equation  becomes 


(3.6) 


In  order  not  to  violate  the  second  law  of  thermodynamics,  the  coefficient  of  the 
last  term  (Tp  ,)  in  equation  3.6  has  to  be  greater  than  0.  By  setting  coefficient  of 
( T y t _ ] ) i 0 , the  nodal  minimum  time  step  yields 


(3.7) 


Node  Type  4 


Figure  3.6  Representation  of  node  type  4 


This  node  is  a capacitance  node  similar  to  node  3 except  it  is  located  at  the 
bottom.(Figure  3.6).  The  temperature  for  the  node  at  the  time  step  ( T1"' ) is 


and  the  maximum  time  step  is 


Figure  3.7  Representation  of  node  S 


Specification 

This  node  is  a capacitance  node  located  at  the  top  center.  It  is  similar  to  node  3 
except  the  left  side  is  insulated  (Figure  3.7).  The  equation  for  the  nodal  temperature  at 


(3.10) 


By  setting  the  coefficient  of 


0,  the  nodal  maximum  time  step  becomes 


(3.1 1) 


Figure  3.8  Representation  of  node  type  6 


Specification 

This  is  similar  to  the  node  type  5 except  it  is  at  the  bottom  (Figure  3.8).  The 
temperature  at  time  step  Tm  is 


(3.12) 


(3.13) 


sfc’«S 
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By  setting  the  coefficient  of  the  last  term  ( T ' ) greater  than  0.  the  maximum  time 


step  becomes 


(3.15) 


Node  Type  8 


Figure  3. 10  Representation  of  node  type  8 


Spec  1 1 t on 

This  is  a capacitance  node  located  at  the  center,  and  its  left  side  is  insulated.  It  is 
surrounded  by  three  capacitance  nodes(Figure  3.10).  The  energy  balance  on  the  node  (i  j) 
is  the  sum  of  the  heat  conducted  from  the  3 adjacent  nodes  is  equal  to  its  internal  energy 
increase.  The  equation  for  the  nodal  temperature  at  time  step  IT'*1)  is 


13.16) 


By  setting  coefficient  of  Tt  . ‘ ' 0,  the  nodal  maximum  time  step  of  this  node  is 


Node  Type  9 


Representation  of  node  type  9 


Specification 

This  is  a capacitance  surface  node  having  a triangular  cross-section  ring  shape, 
with  volume  I v (Figure  3, 1 1).  Heat  convects  from  the  right  boundary  into  the  node 
through  the  sloped  surface  area  (ArlkrllJ),  and  conducts  from  two  capacitance  nodes  on 
the  left  and  top  into  the  node  through  the  average  cross  section  areas  Ar  ^ v andAy 
(see  derivation  in  appendix  A). 


The  energy  I 
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The  nodal  temperature  for  lime  step  ( T p‘ 1 ) be 


(3.20) 


Npde  Type  IQ 


Specification 

This  is  a capacitance  surface  node  having  a shape  of  a triangle  ring  and  the  nodal 
volume  is  It  is  located  at  the  right  top  corner  (Figure  3.12).  Heat  convects  from 
the  right  side  into  the  node  through  the  sloped  area  Ar  and  conducts  from  the  left 
node  through  the  average  cross  section  area.  Node  on  the  top  is  a NCSN  with  convective 
boundary.  The  energy  balance  of  this  node  becomes 


(3.21) 


and  the  equation  for  the  nodal  temperature  for  time  step  ( T*”1 ) is 


(3.22) 


By  setting  coefficient  of  Tltl  . ' ? 0,  the  nodal  maximum  time  step  of  this  node  is 


(3.23) 


Grid  Generation  in  a Finite  Elliptical  Cylinder 
The  equation  of  an  ellipse  is  given  by 


where  ‘2a’  is  the  length  of  the  major  axis,  ‘2b’  is  the  length  of  the  minor  axis.  A two 
dimensional  Unite  difference  model  is  developed  for  a finite  elliptical  cylinder  with  height 


H.  Because  of  ihe  symmetry  of  the  finite  elliptical  cylinder,  only  the  right  half  is  used  to 
develop  the  equations.  Heat  flows  into  the  finite  elliptical  cylinder  from  two  directions: 
along  the  curve  line  drawn  perpendicularly  from  the  node  on  the  ellipse  perimeter  toward 
the  center  (Pham,  1991;  Erdogdu,  1996),  and  along  the  height.  The  equation  of  Ihe  curve 
line  is  given  as 


Figure  3.13  shows  the  generation  of  nodal  network  within  a finite  elliptical 
cylinder.  The  expressions  of  distances  between  nodes,  areas  of  heat  transfer  and  nodal 
volumes  are  presented  in  the  Appendix  B.  More  detail  in  grid  generation  of  an  infinite 
elliptical  cylinder  could  be  obtained  from  Erdogdu(l996). 

Eneruv  Balance  Equations 

From  the  side  view  of  grid  generation  in  a finite  elliptical  cylinder  (Figure  3.13b). 
twelve  different  types  of  nodes  were  generated  within  the  finite  elliptical  cylinder  as 
shown  in  Figure  3. 14.  The  derivation  of  the  energy  balance  equations  and  the  calculated 
maximum  time  steps  for  the  volume  elements  were  types  I to  8 of  the  finite  elliptical 
model  are  the  same  as  those  for  the  cone  frustum  model.  The  following  are  the  developed 
equations  for  the  nodes  types  9 to  1 2 of  the  finite  elliptical  cylinder. 
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Fig  3. 13  Grid  Generation  in  a finite  elliptical  cylinder 


j = 1 tojmax 


Figure  3. 14  Volume  elemem  definitions  in  a finite  elliptical  cylinder 

0 Capacitance  nodes 
O Non-  Capacitance  Surface  Nodes 


Mate  Ins .?  (ELiiK.flliBiiwI  Cylinder) 


Figure  3.15  Representation  of  node  type  9(Finitc  elliptical  cylinder) 


This  is  a capacitance  node  having  a NCSN  with  convective  boundary  on  the  right 
side(Figure3.l5).  The  temperature  in  time  step  t T'”'1 ) is 


(3.26) 


By  setting  coefficient  of  T 


: step  of  this  node  i 


(3.27) 


Node  Tvne  10  (Finite  Elliptical  Cylinder  1 


Figure3-16  Representation  of  node  type  10  (Finite  elliptical  cylinder) 


This  is  a capacitance  node  having  a NCSN  with  convective  boundary  on  its  top 
and  right  sides(Figure  3.16).  The  temperature  in  time  step  ( T'"' ) is 


(3.28) 


By  setting  coefficient  of  T >.  0,  the  maximum  time  step  of  this  node  is 


Node  Type  1 1 (Finite  Elliptical  Cylinder) 


Figure  3 17  Representation  of  node  type  1 1 (Finite  elliptical  cylinder) 


This  is  a capacitance  node,  having  NCSN  nodes  with  convection  boundary 
right  and  bottom  (Figure  3. 1 7),  The  temperature  in  time  step  IT'"1)  is 


By  selling  coefficient  of  T F ? 0,  the  nodal  maximum  time  step  yields 


Node  Tvne  12  (Finite  Elliptical  Cvlinderl 


Figure  3.18  Representation  of  node  type  12  (Finite  elliptical  cylinder) 


Tliis  is  a NCSN,  located  on  the  right  surfacc(Figure  3.18).  Its  temperature  can  be 


Validation  of  the  Developed  Models 


The  2>dimensionaI  finite  difference  model  for  a finite  cylinder,  developed  by  the 
energy  balance  approach  with  non-capacitance  surface  nodes,  was  used  to  validate  the 
cone  and  finite  elliptical  cylinder  models  mentioned  above.  The  accuracy  of  the  cylinder 
model  was  validated  by  comparing  its  results  with  those  of  analytical  solutions  and 
experimental  temperature  measurements  (Erdogdu,  1996).  Details  of  the  development  of 
the  cylindrical  model  are  given  by  Welt  et.al.  (1997). 

Both  the  cone  and  finite  elliptical  cylinder  models  were  assigned  the  same 
dimensions  as  the  cylinder,  except  that  the  large  radius  of  the  cone  could  not  be  set  smaller 
or  equal  to  the  small  radius.  Therefore,  the  large  radius  used  in  the  cone  model  was  set 
slightly  bigger  than  the  small  radius  (Table  3.1)  Grids  and  nodes  in  the  cylinder,  cone  and 


finite  elliptical  cylinder  models  were  generated  so  that  the  temperatures  at  the  same 
locations  could  be  compared. 

The  identical  values  of  thermal  difiusivity  (1.4*10 7 nr  /s),  initial  temperature  (22 
11 C),  heating  temperature  (121  C),  surface  heat  transfer  coefficient  during  heating  (2000 
W/m!  K),  heating  time  (60  min),  cooling  temperature  (20  °C).  surface  heat  transfer 
coefficient  during  cooling  (300  W/  m 7 K)  and  cooling  time  (60  min)  were  used  as  inputs 
to  the  three  models. 

The  center  temperature  histories  and  temperature  distributions  at  different 
locations  simulated  using  the  cone  and  the  finite  elliptical  cylinder  model  were  compared 
with  the  simulations  at  the  some  points  obtained  by  the  validated  cylinder  model. 
Experimental  Validation  with  Acrylic  Materials 

Acrylic  finite  elliptical  cylinders  and  cones  of  different  dimensions  were  used  in  the 
heat  penetration  tests  in  order  to  validate  the  conduction  heat  transfer  models  developed 
The  acrylic  rod  with  a diameter  of  10.24  cm  (Acutech  Plastics.  Reading , PA)  was  formed 
in  to  different  shapes  and  sizes.  The  thermophysical  properties  are  given  in  Table  3.2. 

The  dimensions  of  the  cones  are  shown  in  Table  3.3  and  those  of  the  finite  elliptical 
cylinders  are  shown  in  Table  3,4.  The  acrylic  elliptical  cylinder  # I was  used  in  the  work  of 
Erdogdu  ( 1 996),  and  the  acrylic  elliptical  cylinder  # 2 was  used  in  the  work  of  Kim  and 
Tebteira(l996). 

It  was  important  to  choose  the  locations  where  the  temperatures  were  measured. 
Ideally  the  location  of  the  thermocouple  tip  should  be  exactly  at  the  location  of  the 
numerical  model  node  in  order  to  avoid  interpolation.  After  the  thermocouple  positions 


Table  3.2  Thermophysical  properties  of  the  actylic  rod  used  in  the  experimental  validation 


Properties 


Thermal  Conductivity 
(W/m-K) 


Specific  Heat  (J/kg-K) 
Density  (kg/m*) 


Value 

0.166-0.25 


1464.40 

1170-1200 


e:  Acutech  Plastics,  Reading,  PA. 


Table  3.3  Dimensions  of  the  acrylic  cones  and  locations  of  the  thermocouple  tips  used  in 
the  experimental  validation. 


Cone* I 

Big  Radius  = 5.08  cm. 
Small  Radius  = 0 cm. 
Height  = 12.7cm. 


Conc#2 

Big  Radius  = 5.08  cm. 
Small  Radius  = 2.54  cm. 
Height  = 7.62  cm. 


Cone  #3 

Big  Radius  = 3.81  cm. 
Small  Radius  = 1 .27  cm. 
Height  = 12.7  cm. 


Table  3.4  Dimensions  of  the  acrylic  finite  elliptical  cylinder  and  locations  of  the 
thermocouple  tips  used  in  the  experimental  validation. 


■ assigned,  small  holes  were  drilled  to  insert  i 


touples  at  those 


positions.  The  thermocouples  were  set  using  epoxy  glue  and  sealed  at  the  surface  with 
silicone. 

The  computerized  data  acquisition  system,  Can-Calc  program®,  developed  by  Or. 
M.O.  Balaban  of  the  University  of  Florida  was  used  to  read  temperatures.  Thermocouples 
were  connected  to  the  DAS-TC  A/D  board  from  MetraByte  (Cleveland,  OH).  All 
thermocouples  were  calibrated  with  a mercury-in-glass  thermometer  prior  to  each 
experiment.  In  order  to  obtain  uniform  initial  internal  temperatures,  the  acrylic  shapes 
were  immersed  in  a well  agitated  constant  temperature  water  container  for  at  least  12 
hours.  The  temperatures  at  every  thermocouple  location  within  the  shapes  and  the  water 
temperature  were  recorded  until  the  difference  was  within  0.5°C. 

Heat  penetration  tests  were  carried  out  in  a temperature  controlled  water-bath 
with  adjustable  water  agitator  set  at  the  highest  speed  (Fischer  Scientific,  Pittsburgh,  PA). 
Once  a constant  water  temperature  was  reached,  the  shape  was  promptly  transferred  into 
the  water  bath.  The  water  temperature  and  temperatures  within  the  acrylic  model  were 
recorded  every  10  seconds.  In  experiments  where  variable  water  temperature  was 
needed,  ice  water  was  used  to  adjust  the  water  temperature.  Both  constant  and  variable 


The  water  temperatures  used  in  the  simulations  were  the  actual  temperatures 
recorded  during  the  heat  penetration  tests.  For  each  shape,  optimum  number  of  nodes  wt 


found  by  repeating  calculations  with  increasing  the  number  of  nodes  until  changes  i 


calculated  temperatures  were  negligible.  The  predicted  temperature  histories  at  every 
thermocouple  location  were  compared  with  experimental  values. 

Calculation  for  Retention  of  Quality  Factors 

The  calculation  for  retention  of  any  quality  factors  makes  use  of  the  temperature 
predicted  at  every  volume  element  after  each  iterations  during  the  numerical  solution  of 
the  conduction  heat  transfer  equations  Since  the  D value  of  a quality  factors  changes  as  a 
function  of  temperature  and  Z value,  it  was  updated  for  every  volume  element  as 

Dr  - lO113^  <333> 

Using  the  updated  D value  (DT),  the  retention  after  a time  change  of  At,  could  be 
calculated  os 


= N,  l0'  41,0,1  (3.34) 

An  initial  retention  of  100%  was  used  and  the  percent  volume  average  retention  after 
every  time  interval  was  calculated  by  volume-averaging  the  retention  throughout  the 

Calculation  of  Sterilizing  Value  (F_) 

The  F0  is  the  integrated  lethality  calculated  by  the  general  method  which  is  based 
on  the  actual  temperature  at  the  cold  point  for  the  entire  process 


where  T is  temperature  at  the  cold  point  at  time  t,  and  the  Z value  was  assumed  to  be 
I0"C. 

The  position  of  cold  point  for  the  finite  elliptical  cylinder  and  cylinders  are  at  the 
geometrical  center.  In  a cone  frustum,  the  position  of  the  cold  point  depends  on  the  shape 
and  size.  Therefore,  it  is  necessary  to  search  for  a cold  point  position  and  calculate  its  F„. 
The  cold  point  was  searched  for  by  comparing  temperatures  of  nodes  along  the  center  line 
in  the  axial  direction.  A fine  grid  was  required  to  obtain  an  accurate  cold  point  position. 


The  objective  (unction  was  the  volume  average  quality  retention,  expressed  as 


The  constraint  of  the  optimization  procedure  was  a target  F„  value  at  the  cold 
point  including  both  heating  and  cooling  period  (t  h , t ,) 


Optimization  of  Quality  Retention 


(3.36) 


(3.37) 


- a.  ; » - a, 


A !\  value  of  3 minutes  was  used.  Although  this  value  is  less  than  those  generally 
used  in  the  industry,  it  shortened  that  calculation  times  in  this  study 

Quality  factors  with  different  combinations  of  D and  Z values  were  used  in  this 
study  Influential  variables  including  thermal  diflusivity  of  food  product,  surface  heat 
transfer  coefficient,  as  well  as  shapes  and  sizes  of  the  cones  and  Unite  elliptical  cylinder 
were  used  in  the  investigation  to  study  their  effect  on  those  factors  The  values  chosen  for 
these  variable  are  summarized  in  Table  3.5, 

The  food  products  with  uniform  initial  temperature  of  22 " C were  assumed  to  be 
processed  at  different  constant  retort  temperatures  ranging  from  1 10  to  140  “C  at  2 “C 
intervals.  Afler  heating  at  the  specified  temperature,  the  food  was  then  cooled  at  a 
constant  cooling  temperature  of  20  ° C until  the  center  temperature  was  less  than  or  equal 
to  50-C. 

The  contribution  of  F„  during  cooling  should  not  be  neglected  for  a conduction 
heating  product.  Depending  on  the  cold  point  temperature  at  the  start  of  cooling,  size  of 
the  container  and  other  factors  in  a conduction  heating  situation,  the  contribution  of 
cooling  to  F„  could  be  as  much  as  40%  of  the  total  (Datta  et  al.,1 986).  In  this  study,  the 
computer  program  developed  for  simulation  of  quality  retention  can  determine  the  exact 


Table  3.5  Variables  used  in  the  quality  retention  studies 


Table  3.6  Dimension  ofthe  cones  used  in  the  quality  retention  studies 


Dimension 
Big  Radius  (cm.) 
Small  Radius! cm.) 


Volume  (cm*) 
Volume  / Surface 


Cone  U I 
4.50 
2.75 


Cone  n 


5.00 


273.90  326.30 

1.23  1.03 


Cone  #3 

4.00 

2.00 
11.50 

337.20 

282.80 

1.19 


Table  3.7  Dimension  of  the  finite  elliptical  cylinder  used  in  the  quality  retention  studies 


Semi  Major  Axis,  a (cm.) 


Semi  Minor  Axis,  b (cm.) 


Volume  / Surface 


308.71 

321.01 


healing  and  cooling  times  al  the  given  heating  temperature  to  obtain  the  specified  F„  value 
and  calculate  nutrient  retention  at  the  end  of  that  process. 

The  exact  time  to  stop  heating  and  start  cooling  in  order  to  obtain  the  desired  F„ 
value  was  found  as  follows:  At  the  end  of  a given  time  step  during  heating,  the  F0  for 
heating  was  calculated,  then  cooling  was  simulated  until  the  cold  point  reached  50  "C.  and 
the  contribution  of  cooling  to  F0  was  calculated.  If  the  total  F„  (heating  + cooling)  was  less 
than  the  desired  F„,  heating  could  still  proceed.  Otherwise,  a bisection  method  was  used 
within  the  last  time  step  to  find  the  exact  time  such  that  the  total  F„  would  not  be  larger 
than  0,0001  min  compared  to  the  desired  F„ 

The  percent  quality  retention  for  each  processing  condition  was  plotted  against 
processing  temperatures.  The  optimum  temperature  giving  the  maximum  quality  retention 
svas  located  by  a spline  interpolation. 


CHAPTER  4 

RESULT  AND  DISCUSSION 


Heal  Transfer  Model  Validation 


The  center  temperature  histories  of  a cylinder  predicted  by  the  cone  model  and 
cylinder  model  were  compared,  and  are  shown  in  Figure  4, 1 . The  agreement  was  well 
within  a variation  of  0.6  °C  lor  the  entire  process.  Because  of  the  limitation  of  the  cone 
model,  the  big  radius  used  in  the  cone  model  was  slightly  bigger  than  the  radius  of  the 
cylinder  used  in  the  cylinder  model.  This  difference  in  dimensions  caused  the  center 
temperature  predicted  using  the  cone  model  to  be  slightly  lower  during  heating  and  higher 
during  cooling  compared  to  that  predicted  by  the  cylinder  model. 

Afler  heating  and  cooling,  the  temperature  distributions  of  the  right  bottom  quarter 
of  the  cylinder  obtained  from  the  cone  model  and  the  cylinder  model  were  compared  in 
Table  4.1.  The  results  show  very  good  agreement  of  temperatures  for  all  locations  with 
variations  within  0.6 "C  for  all  locations  indicated  that  the  cone  model  could  predict 
temperature  distribution  of  a cylinder-like  object  as  accurately  as  the  validated  cylinder 


Figure  4. 1 Comparison  of  center  temperature  of  307*409  can  predicted  using  the  cone 
model  and  cylinder  model 


However,  using  (he  cone  model  for  predicting  temperature  of  a cone  which  has  a 
shape  close  to  a cylinder  is  not  recommended.  In  the  grid  generation  of  the  cone  model, 
nodes  in  the  cone  height  were  automatically  set  equally  to  the  difference  between  nodes  in 
a big  radius  and  a small  radius.  If  this  difference  is  very  small,  the  elements  in  the  sloped 
area  become  very  small.  Since  the  forward  difference  explicit  procedure  was  used  for  the 
finite  difference  approximations,  the  maximum  time  step  became  smaller  when  the  grid 
size  is  decreased  resulting  in  longer  computer  calculation  times 
Mathematical  Validation  pf  the  nnijeEHiplical  Cylinder  Model  with  Existing  Cylindrical 
Model 

To  determine  the  accuracy  of  the  developed  model,  the  center  temperature 
histories  of  a cylinder  (Eccentricity  = 0)  predicted  by  the  finite  elliptical  cylinder  model 
were  compared  with  that  obtained  from  the  validated  cylinder  model  for  the  same 
conditions  of  size,  thermal  diffusivities,  initial  temperatures,  outside  temperatures  and 
surface  heat  transfer  coefficients.  The  predicted  temperatures  from  the  two  models  were 
identical  throughout  the  heating  and  cooling  periods  (Figure  4.2). 

The  final  temperature  distributions  obtained  from  the  two  models  were  also 
identical  (Table  4.2).  This  validation  illustrated  that  there  was  no  difference  in  the 
temperature  prediction  for  a finite  cylinder  using  the  finite  elliptical  cylinder  model  and  the 
cylinder  model  However,  at  the  same  accuracy,  the  finite  elliptical  cylinder  model  was 
slower  since  the  elements  created  in  the  finite  elliptical  cylinder  model  were  smaller  than 
those  in  the  cylinder  model. 


Figure  4.2  Comparison  of  center  temperature  of  307*409  can  predicted  using  the  finite 
elliptical  cylinder  model  and  the  cylinder  model 


75 

Experimental  Validation  with  Acrylic  Materials 

An  acrylic  material  was  used  as  a model  for  experimental  validation  of  the 
developed  heat  transfer  models  because  it  is  homogeneous;  heats  by  conduction;  has 
known  thermal  properties  that  could  be  obtained  from  the  manufacturer,  and  is  easily  to  be 
shaped  into  different  geometries.  The  experimental  temperatures  at  the  different  locations 
within  2 finite  elliptical  cylinders  and  the  3 acrylic  cones  were  compared  with  the 
simulation  temperatures  (Figures  4.3-  4.12).  The  corresponding  heat  transfer  coefficients 
and  thermal  properties  used  in  the  experiments  are  included  in  those  figures.  The 
dimensions  and  the  thermocouple  locations  arc  given  in  Table  3.3.  The  experimental  data 
can  be  seen  in  Appendix  C (Table  C.  I-C.  10). 

The  experimental  results  agree  well  with  the  simulations  for  all  locations  of  every 
geometrical  shape  in  response  to  either  constant  or  variable  outside  temperatures  The 
heat  transfer  coefficient  was  selected  by  a search  procedure  to  obtain  a minimum  absolute 
difference  between  simulated  and  experimental  temperatures  for  all  thermocouple 
locations.  In  the  search  procedure,  the  values  of  thermal  conductivity  and  diffusivity  (1.9 
W/m  K and  1 .09  mVscc  respectively)  values  given  by  the  manufacturer  were  used  for  all 
conically  shaped  acrylic  model.  Surface  heat  transfer  coefficient  values  used  to  obtain 
good  agreement  between  experiment  and  simulation  for  all  points  were  within  75-150  W/ 
m'  K.  These  values  were  in  the  same  range  of  values  reported  by  Lebowitz  and  Bhowmik 
( 1 989)  who  experimentally  determined  h value  of  1 86  ± 54  W/  m1  K for  thermal 
processing  of  thin  retortable  pouches  under  120"  C pressurized -water  with  forced 
convection.  In  a study  of  heat  transfer  in  a bowl  which  has  a similar  shape  to  a cone. 
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Sheen  et.al.  (1993)  used  an  h value  of  700  W / m'  K to  simulate  temperature  within  a 
bowl  (big  radius  4.2  cm,  small  radius  2.8  cm  and  height  S.7  cm)  heated  in  a circulating 
superheated  water-type  retort  at  constant  temperature  of  I20°C. 

The  heat  transfer  coefficient  required  for  a good  agreement  of  the  finite  elliptical 
cylinder#  1 was  1200-1300  W/  m: " C'  and  the  finite  elliptical  cylinder#2  was  approximately 
400  W m’ 0 C.  According  to  Singh  and  Heldman  ( 1 993),  the  approximate  value  of 
convective  heat  transfer  coefficient  value  of  forced  convection  water  (not  boiling  )cou!d 
vary  from  SO- 1 0,000  W/  nr  " C Chau  and  Snyder  (1988)  measured  the  h value  of  a 
shrimp  shaped  aluminum  casting  (similar  to  a cone  shape),  heated  in  agitated  hot  water. 
They  reported  the  h values  decreased  from  5472  to  3478  W/  m:  K when  the  radius 
increased  from  6,35  to  8.75  mm.  This  study  should  support  the  contention  that  the 
variation  of  the  heat  transfer  coefficient  in  this  study  is  mostly  due  to  differences  in  shape 
and  size  of  the  acrylic  models. 

It  could  be  concluded  that  the  close  agreement  of  replicate  experimental  and 
simulated  temperature  profiles  of  3-4  points  within  different  shapes  and  sizes  of  the  cone 
and  finite  elliptical  cylinder  models  provided  sufficient  evidence  to  validate  their  accuracy. 

S.tu<!y_9.f  Op.tiinizaii.ojuC<>ality  Kerennon 

The  cone  and  finite  elliptical  cylinder  models  were  further  modified  to  calculate  the 
volume  average  quality  retention  of  foods  during  thermal  processing.  A number  of 
simulations  were  performed  with  primary  emphasis  on  a dear  understanding  of  the 
influence  of  various  parameters  on  the  maximum  quality  retention.  Six  Z values  (10, 20, 


30, 40,  50  and  60  "C  ) and  six  D„,  values  (40, 80, 120,  160, 200  and  240  minutes), 
resulting  in  36  combinations  of  quality  factors  covering  the  broad  range  of  kinetic 
parameters  found  in  the  literature  were  considered.  The  rate  of  thermal  destruction  of 
these  quality  factors  was  assumed  to  follow  first  order  kinetics.  The  effects  of  processing 
conditions  (retort  temperature,  surface  heal  transfer  coefficient,  and  container 
geometries ) and  the  thermal  diffusivity  of  foods,  as  summarized  in  Table  3.5,  was 
systematically  investigated.  Three  equal-volume  conically  shaped  containers  with 
different  geometries  and  two  typical  commercial  finite  elliptically  shaped  containers  with 
dimensions  given  in  Table  3.6  and  3.7  were  used  in  the  simulations.  Internal  heat  transfer 
was  conduction  only.  It  is  assumed  that  the  temperature  at  the  surface  of  the  object  being 
heated  is  constant  and  thermal  difliisivity  of  the  object  is  constant  over  the  temperature 
ranges  under  consideration. 

In  this  study,  the  cold  point  at  any  given  time  was  defined  as  the  point  of  the  least 
accumulated  F„  at  that  time.  Because  conduction  heating  food  in  a conically  shaped 
container  is  heated  asymmetrically,  the  position  of  the  cold  point  depends  on  the  container 
geometry. 

Figure  4.13  shows  the  simulation  of  temperature  distribution  along  the  axis  of  the 
3 conically  shaped  containers  during  heating  for  30  minutes  at  10  minute  intervals  in  a 
constant  1 20"C  heating  medium  under  an  infinite  heat  transfer  coefficient  condition.  The 
position  index  indicates  the  location  along  the  axis  from  the  small  radius  (position  index 
= 0.0  )to  the  big  radius  (position  index  =1.0).  The  effect  of  food  thermal  diffusivity  at 


1.2, 1.4  and  1 .6  nr/sec  on  the  cold  point  location  is  also  shown  in  the  same  figure.  The 
location  of  the  cold  point  did  not  change  over  time  during  the  heating  of  1 0 to  30  minutes 
regardless  of  cone  geometry  or  thermal  diffusivity.  The  position  indices  of  the  cold  point 
for  the  cone  #1,  #2  and  S3  were  0.58, 0 50and  0.71,  respectively 

The  cold  point  location  during  heating  for  all  conical  shapes  did  not  change  when 
the  heat  transfer  coefficient  changed  from  300  W m1  /K  to  infinite  h (Figure  4. 14),  With 
an  infinite  heat  transfer  coefficient  condition,  the  temperatures  at  the  position  indices  of 
0.0  and  1 .0  were  equal  to  the  heating  medium  temperature.  As  the  heat  transfer 
coefficient  decreased,  the  temperature  of  all  points  along  the  rotation  axis  decreased. 

The  simulation  results  indicated  that  the  position  of  the  least  F„  may  change  alier 
cooling  has  started  depending  on  cone  geometry  and  may  not  be  at  the  cold  point  during 
heating.  Since  the  cold  point  is  cooled  more  slowly,  the  contribution  of  F„  is  still  at  a high 
rate  and  the  total  accumulated  F„  at  this  point  may  be  higher  than  the  total  F„  of  other 
points  along  the  axis.  In  order  to  find  the  processing  time  that  gives  the  desired  F„  value,  it 
is  important  to  keep  track  of  the  position  of  the  least  F„  for  the  entire  process  and 
compare  it  with  the  target  F„  to  ensure  the  safety  of  the  product. 

Processing  Time  for  Equivalent  Lethality 

The  simulations  assumed  that  food  was  heated  from  the  initial  temperature 
( 22  °C)  and  cooled  down  until  the  cold  point  temperature  was  less  than  50  °C.  The 
lethality  (F0)  of  3 minutes  which  could  ensure  safety  from  Clostridium  bolulinum  was  used 
as  a constraint  in  the  simulations. 


The  exact  time  for  heating  and  cooling  in  order  to  obtain  the  desired  F„  value  was 
found  such  that  the  total  F„  would  be  not  exceed  the  desired  F„  by  more  than  0.0001  min 
The  total  processing  time  was  the  sum  of  the  heating  and  cooling  time.  The  complete 
simulated  heating,  cooling  and  total  processing  time  responding  to  different  processing 
temperatures  for  both  conical  and  elliptical  containers  were  summarized  in  Appendix  D 
(Table  D- 1 and  D-2).  The  total  time  and  corresponding  retort  temperature  required  for 
processing  different  container  geometries  were  compared  in  Figure  4. 1 5.  Any  point  on 
these  curves  represents  a thermal  process  that  produces  the  same  target  F„  of  3 minutes 
with  the  error  less  than  1%.  Since  the  destruction  rate  of  microorganisms  has  an 
exponential  relationship  with  temperature,  as  expected,  the  relationship  of  retort 
temperature  and  processing  time  is  also  an  exponential  function.  When  a low  retort 
temperature  was  used  especially  when  it  was  for  below  the  reference  temperature  (121.1 
C),  the  longer  processing  time  was  required  to  reach  a target  F„  because  a very  small 
lethality  is  contributed  at  each  time  interval.  As  the  processing  temperature  was  raised 
closer  to  the  reference  temperature,  the  processing  time  was  significantly  shortened  and 
then  it  remained  almost  constant. 

Among  the  conical  containers  of  equivalent  volume  (approximately  336  cm’),  the 
cone  # 2 which  had  the  smallest  volume-lo-surfacc  ratio  (1 .03  ) required  the  shortest 
processing  time,  followed  by  the  cone  #3  and  ff  1 with  ratios  of  1 . 19  and  1 .23, 
respectively. 

Comparing  between  the  finite  elliptical  cylinder#  1 and  #2,  having  volume  of  5 1 5 
and  308  cm',  as  expected,  the  finite  elliptical  cylinder  # 1 required  longer  processing  time. 


It  is  interesting  to  note  that  the  finite  elliptical  cylinder  H I , which  has  a bigger  volume  and 
equal  volume/surfacc  ratio(l  19)  compared  to  the  cone  H3,  required  less  processing  time. 
This  is  because  the  finite  elliptical  cylinder  #1  has  a large  symmetric  top  and  bottom 
surface  and  a thin  profile  that  makes  heat  transfer  to  the  cold  point  faster.  The  simulation 
results  lead  to  the  conclusion  that  volume,  surface  area  or  even  volume/surface  area  may 
not  be  enough  for  comparing  rates  of  heat  transfer  in  irregular  containers.  Simulation 
should  be  performed  case  by  case  in  order  to  make  a comparison. 

Effects  of  thermal  diffusivity  (1.2,  I.dand  1 .6  si  O'1  nr/s)  and  surface  heat  transfer 
coefficient  (h-  300  W mV  K and  infinite  h ) on  total  processing  time  of  foods  packed  in 
the  cone  #1  is  illustrated  in  Figure  4. 16.  As  supported  by  Figure  4. 13a.  increasing 
thermal  diffusivity  results  in  (aster  heating  at  the  cold  point,  which  then  accelerates  the 
temperature  uniformity  throughout  the  container  and  shortens  the  processing  time. 
Increasing  the  heat  transfer  coefficient  also  increases  the  rale  of  heat  transfer  by  increasing 
the  surface  temperature.  In  the  case  of  cone  HI,  foods  processed  under  infinite  heat 
transfer  coefficient  conditions  require  1 3-24  % less  processing  time  than  those  processed 
under  the  finite  heat  transfer  coefficient  condition. 

It  should  be  noted  that  the  contribution  of  F„  during  the  cooling  of  conduction 
heating  food  is  important  especially  when  food  is  processed  at  high  retort  temperatures.  In 
such  cases,  the  temperature  at  the  cold  point  is  high  when  cooling  starts  and  therefore  it 
takes  longer  time  to  cool  down.  The  ratio  of  heating  to  cooling  time  was  almost  1:1 


Processing  Tertperalure  ( Q 


Figure  4. 1 S Simulated  retort  temperatures  and  total  processing  times  of  conduction  heated  foods 
having  alpha  = 1 .4  nr/s  packed  in  different  geometric  containers  giving  equal  F„  value  = 3 minutes 
( based  on  z value  = 1 0 minutes.  Initial  temperature  = 22  °C  final  temperature  = 50  "C,  infinite 
heat  transfer  coefficient  condition) 
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Figure  4. 16  Simulated  retort  temperatures  and  total  processing  times  of  conduction  heated 
foods  having  different  thermal  diffuisivity  packed  in  the  conical  shaped  container  "1  processed 
under  finite  and  infinite  heat  transfer  coefficient  conditions,  giving  equal  F.,  value  - 3 minutes 
( based  on  z value  - 1 0 minutes,  Initial  temperature  = 22  "C  final  temperature  - 50  ”C) 


Investiualion  of  Factors  Affectinu  Oualiiv  Retention  and  Oniinmm  Processing 


All  results  of  simulated  overall  retention  of  quality  factors  with  different  D„,  and  Z 
values,  in  conduction  heating  foods  with  thermal  diffusivity  1.2,  1.4  and  1 .6  nr/s  packed 
in  2 different  conical  shaped  containers  and  2 different  elliptical  shaped  containers 
processed  with  constant  temperature  from  1 10  to  140  C using  finite  and  infinite  surface 
heat  transfer  coefficients  are  presented  in  Appendix  E (Tables  EI-E30).  Optimum 
processing  temperatures  for  maximizing  quality  retention  of  different  geometries,  thermal 
diffusivity  and  heat  transfer  coefficients  are  summarized  in  Tables  4. 3-4.7. 

Z Value 

The  effect  of  changing  the  Z value  of  quality  factors  across  the  range  of  10-60  “ C 
at  1 0 "C  intervals  and  processing  temperatures  (110- 1 40  ° C)  on  the  volume  integrated 
retentions  is  illustrated  in  Figure  4. 1 7.  These  simulations  were  performed  with  foods  of 
thermal  diffusivity  1 .4  x 10  7 nr/s  packed  in  the  cone  # I and  the  D,.,  value  of  the  quality 
faciors  was  fixed  at  200  minutes.  The  curve  for  depicting  the  effect  of  increasing 
temperature  at  each  Z value  had  one  optimum  temperaIure(T  <)  resulting  in  a maximum 
quality  retention.  The  quality  factor  with  a lower  Z value  had  a lower  T^.  The  results  also 
indicated  that  percent  retention  of  the  quality  factors  with  a low  Z value  are  highly 
sensitive  to  a departure  from  the  optimum  temperature.  Since  these  quality  factors  are 
very  sensitive  to  temperature  changes,  their  destruction  rales  changed  dramatically  when 
the  temperature  changed.  When  the  processing  temperature  is  higher  than  the  optimum 
temperature,  the  effect  of  high  temperature  destruction  overcomes  the  benefit  of  the  high 
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temperature  short  time  process  and  the  integrated  volume  retention  drops  sharply.  From 
Figure  4. 1 7,  at  a Z value  = 10  “C,  the  optimum  temperature  giving  the  maximum 
retention  of  94.4%  occurs  at  1 10  "C.  The  retention  drops  from  94.4%  to  37.1  % upon 
increasing  the  retort  temperature  from  1 10  to  140  °C. 

When  the  Z value  increases,  the  reduction  rale  of  quality  factors  is  more  insensitive 
to  changes  in  temperature.  The  optimum  temperature  for  maximizing  the  nutrients  moved 
toward  higher  temperatures.  However,  as  the  Z value  increases,  the  profile  of  the  curve 
flattens  and  the  percent  quality  retention  becomes  less  sensitive  to  deviations  from  the 
optimum  temperature.  For  these  simulations,  the  quality  factor  ofZ  = 60  °C  had  an 
optimum  temperature  at  132.4'C  (Table  4.3  b)  and  the  percent  retention  at  this  point  was 
6I.6S  %.  The  results  from  Table  E-5  showed  that  using  processing  temperatures  ranging 
from  1 22-1 40  “C,  the  corresponding  quality  retentions  were  not  significantly  different  ( % 
retention  = 59,9-61 .7).  The  difference  from  the  maximum  retention  was  only  ± 1%.  The 
attempt  to  find  the  optimal  temperature  to  maximize  these  quality  factors  may  not  be 
practical  since  the  improvement  is  minor. 

-EUi- 

The  effect  of  changing  the  D„r  value  on  volume  average  quality  retention  was 
investigated.  A plot  of  the  percent  quality  retention  as  a function  of  processing 
temperature  and  DIrt  value  ranging  from  40  to  240  minutes  at  40  minute  intervals,  with  a Z 
value  constant  at  30  °C , is  presented  in  Figure  4. 18.  These  simulations  were  performed 
assuming  food  with  thermal  diffusivity  of  1.4x10  '7  nr/s  packed  in  the  cone  #1.  The 
simulation  results  confirm  the  finding  of  Teixeira  et.al.(1969)  and  Silva  el.al.(!992)  for 


Figure  4.17  Effect  of  Z value  (°C)  and  processing  temperature  (°C)  on  estimated 
retention  of  quality  factor  with  D„i  = 200  min  in  conduction  heated  food  having 
thermal  diffusivity  1.4  x I O'7  m2/s  packed  in  the  conically  shaped  container  ti ! . 
processed  under  infinite  heat  transfer  coefficient  condition 


Figure  4.18  HfTecl  of  D„r  value  and  processing  temperature  (°C)  on  estimated 
retention  of  quality  factor  with  Z value  = 30  °C  in  conduction  heated  food  having 
thermal  diffiisivity  1.4  x I O'7  m2/s  packed  in  the  conically  shaped  container#!, 
processed  under  infinite  heat  transfer  coefficient  condition 


conduction  heated  food  in  cylindrical  container,  quality  retention  increases  when  D„,  value 
increases  at  all  processing  temperatures.  From  Table  4.3  b,  the  maximum  quality  retention 
increases  from  24.45  to  78.00  % when  Dld  increases  from  40  to  240  minutes.  The  D„, 


minutes,  the  corresponding  optimum  temperatures  were  almost  constant  at  approximately 
116  °C.  However,  the  quality  factor  with  low  Devalues  hod  a slightly  lower  optimal 
temperature.  The  maximum  percent  retention  at  D^  = 40  and  80  minutes  occurred  at  the 
temperatures  of  1 19.95  and  1 1 7.42  "C.  respectively.  These  results  agree  with  Silva  et.al. 
(1992),  who  also  found  that  the  Dnl  value  lower  than  150  minutes  had  an  efTect  on  the 
optimal  temperature  of  a conduction  heating  food  in  a cylindrical  can. 

With  a low  Dlt(  value  and  a high  Z value,  the  temperature  required  to  maximize 
quality  retention  is  high.  However,  retention  of  quality  factors  having  such  kinetic 
parameters  is  extremely  limited  at  any  temperatures..  The  maximum  retention  for  that  of 
0^=40  minutes  and  Z = 60  °C  was  only  10.83%.  Similar  trends  regarding  optimum 
temperature  as  a function  of  were  observed  for  different  values  ofZ,  thermal 
diflusivity,  and  other  processing  conditions. 


The  effect  of  thermal  diflusivity  of  1.2,  I 4 and  1.6  xIO'’  mVsec,  which  is  the  range  for 
most  foods,  on  the  volume  integrated  quality  retention  was  investigated.  The  3-D  plot  of 
percent  retention  as  a function  of  thermal  diflusivity  and  processing  temperature  is 
presented  in  Figure  4. 19.  The  simulations  were  performed  with  foods  packed  in  the  cone 
H I under  an  infinite  heat  transfer  coefficient  for  a quality  factor  with  Ddl  = 200  minutes 


value  had  little  influence  on  opi 


. When  D,d  ranges  from  120  to  240 


Figure  4 .19  Effect  of  thermal  diffusivity  (x  I O'7  mVs)  and  processing  temperature  (°C) 
on  estimated  retention  of  quality  factor  with  Z value  = 30  °C  and  D,-(  value  = 200  min 
packed  in  the  conically  shaped  container  If  I,  processed  under  infinite  heal  transfer 
coefficient  condition 


and  Z = 30  °C,  which  corresponds  to  a range  of  kinetic  parameters  for  vitamins  (Table 
2.2).  Foods  with  a higher  thermal  diflusivity  always  retain  higher  quality  for  any 
processing  temperature  because  they  heat  faster  resulting  in  shorter  process  time.  The 
maximum  percent  retentions  increased  from  67.47  to  70.08  and  then  to  72.31  as  thermal 
diflusivitics  increased  from  1 2 tol.4  tol.6,  respectively.  Although  the  thermal 
diffusivity  had  little  influence  on  the  resulting  optimal  temperature,  the  increase  in  optimal 
temperature  when  thermal  diffusivity  increases  indicates  that  the  faster  heating  rate  of 
foods  with  high  thermal  diflusivity  allows  them  to  respond  better  for  the  benefit  of  the 
HTST  processing  to  retain  their  qualities. 

This  increase  in  quality  factor  retention  for  food  with  higher  thermal  diflusivity  was 
also  true  for  all  kinetic  parameters  regardless  of  container  shapes  and  sizes,  and 
processing  temperatures.  This  is  because  the  quality  improvement  is  a benefit  of  the  fester 
heat  transfer  rate  in  foods  with  higher  thermal  diflusivity  which  allowed  the  temperature 
distribution  of  the  food  to  be  uniform  faster  and  reduce  processing  time  for  all  processing 
temperature. 


The  infinite  heat  transfer  coefficient  assumption  is  valid  when  (bod  is  packed  in  a 
metal  container  and  heated  with  condensing  steam.  When  food  is  packed  in  a non-metal 
material  and/or  heated  with  water  or  steam/air  mixtures,  a finite  surface  heat  transfer 
coefficient  is  considered.  In  this  study,  h ■ 300  W/  nr  K was  used  for  finite  h;  the 
volume  average  quality  retention  was  compared  with  those  using  the  infinite  h condition. 


The  surface  temperature  of  foods  heated  under  the  infinite  surface  heat  transfer 
coefficient  condition  will  be  immediately  raised  to  the  medium  temperature  (Figure  4.14). 
High  surface  temperature  is  a driving  force  tor  increasing  the  rate  of  heat  transfer 
throughout  the  container,  resulting  in  a reduction  of  processing  time.  However,  the 
simulation  result  indicated  that  the  volume  integrated  quality  retention  was  insensitive  to 
h value.  Changing  h from  finite  to  infinite  changed  less  than  ± 3%  of  percent  retention 
regardless  of  shape  of  a container,  thermal  diffusivity  of  the  food  and  /or  the  kinetic  data 
of  food  quality  factors  (Appendix  E,  Table  E 1-E  30). 

Even  though  the  effect  of  h is  very  small,  the  simulation  result  indicated  that  the 
quality  factors  with  low  Z value,  processed  at  high  temperature  with  infinite  heat  transfer 
coefficient,  had  lower  percent  retention  than  those  processed  with  finite  h (Figure  4.20). 
This  is  because  with  infinite  h,  the  food  surface  temperature  is  high  therefore  the  quality  at 
the  surface  is  mostly  destroyed  by  heat  before  the  target  lethality  is  reached  at  the  cold 
point.  Since  the  volume  at  the  surface  constitutes  a large  proportion  of  a container,  as  the 
overall  volume  retention  was  integrated,  a lower  retention  was  obtained. 

Quality  factors  with  high  Z values  had  a higher  retention  with  an  infinite  h than  a 
finite  h.  This  is  because  their  thermal  destruction  rate  does  not  change  much  with  higher 
temperature,  and  they  benefit  from  the  shorter  processing  time  using  infinite  heat  transfer 
coefficicnt..lt  was  also  found  that  the  optimum  temperatures  for  foods  using  finite  heat 
transfer  coefficient  were  always  slightly  higher  than  those  for  infinite  heal  transfer  cases. 
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Figure  4.20  Effect  of  heat  transfer  coefficient  on  simulated  retention  of  quality  factors 
with  Devalue = 200  minutes  in  conduction  heated  food  having  thermal 
diffusivity  1.4x  lO^mVs  packed  in  the  conically  shaped  container  41 


Container  Geometry 


The  effect  of  different  geometries  of  conical  and  elliptical  containers  on  volume 
integrated  quality  retention  is  presented  in  Figure  4.2 1 . The  simulations  were  performed 
with  foods  with  thermal  diflusivity  of  1 .4  x I O'1  nr'/s  using  infinite  heat  transfer  coefficient, 
for  a quality  factor  with  D„f  = 200  minutes  and  Z = 30  °C.  The  level  of  quality  retention 
could  be  related  to  a total  processing  time  required  for  each  container  to  achieve  equal 
lethality  (Figure  4.15).  The  results  indicate  that  as  foods  are  packed  in  a container 
requiring  less  processing  time,  they  retain  higher  quality  at  all  processing  temperature. 
This  applies  for  all  quality  factors,  in  any  processing  condition  within  the  scope  of  this 
study  Although  the  containers  had  very  close  optimum  temperaturcs(Figure  4.2 1 ),  the 
containers  that  require  shorter  processing  time  trend  to  have  higher  optimum 
temperatures. 

Quality  reduction  in  foods  during  thermal  processing  is  unavoidable.  However, 
with  a proper  container,  the  quality  loss  could  be  reduced  while  ensuring  product  safety 
The  containers  that  allow  rapid  heat  transfer  and  temperature  uniformity  lead  to  a shorter 
processing  time.  As  a result,  foods  retain  more  quality/nutritional  value,  while  using  less 
energy.  The  developed  heat  transfer  models  are  useful  tools  for  evaluating  quality 
retention,  assisting  the  food  processor  in  choosing  the  appropriate  container  geometiy  and 
optimum  processing  temperature  to  achieve  a better  product  quality. 


CHAPTER  5 
CONCLUSION 


i The  conduction  heat  transfer  models,  developed  with  the  finite  difference 
method  based  on  energy  balance  approach  for  any  given  shape  and  size  of  a conical  and  a 
finite  elliptical  shaped  container,  can  accurately  predict  the  temperature  distribution  in 
response  to  either  constant  or  variable  surface  temperature. 

2.  The  developed  models  can  be  further  modified  to  predict  the  volume  integrated 
retention  of  any  quality  factor  following  first  order  reaction  kinetics,  under  the  constraint 
of  achieving  a given  target  lethality. 

3.  Quality  retention  for  factors  with  low  Z values  are  more  sensitive  to  changes  in 
process  temperature  away  from  the  optimum  temperature  than  those  with  high  Z values. 
The  exact  determination  of  the  optimum  temperature  for  the  quality  factors  with  Z less 
than  50  "C  may  not  be  necessary  since  only  a small  improvement  in  percent  retention  will 
be  obtained. 

4.  Quality  retention  increases  when  D..,  value  increases  at  all  processing 


temperatures.  Although  li,,.  value  had  little  influence  on  optimum  temperatures,  the 
quality  factor  with  low  D„,  values  with  high  Z value  had  a slightly  lower  optimal 


5.  Foods  wiih  a higher  thermal  difTusivity  require  a shorter  processing  time  and 
retain  higher  quality  for  any  processing  temperature.  The  optimal  temperature  slightly 
increases  when  thermal  diffusivity  increases  from  1 .2  to  1 .6  x 1 0 m’/s. 

6.  Even  though  foods  processed  under  infinite  heat  transfer  coefficient  conditions 
require  a shorter  processing  time  than  those  processed  under  finite  heat  transfer  coefficient 
conditions,  the  level  of  quality  retention  after  processing  under  these  conditions  is  hardly 
affected.  Surface  heat  transfer  coefficient  had  a small  influence  on  the  optimum  processing 
temperature. 

7.  A container  that  allows  rapid  heat  transfer  results  in  foods  retaining  more 
qualuy/nuiritional  value  at  all  processing  temperature.  Although  the  containers  withing 
the  scope  of  this  study  had  very  close  optimum  temperatures,  the  containers  that  require 
shorter  processing  time  trend  to  have  higher  optimum  temperatures  Heat  transfer  models 
provide  the  food  processor  an  effective  tool  in  choosing  the  appropriate  container 
geometry  to  achieve  a better  product  quality. 


CALCULATION  OF  NODAL , 


; FRUSTUM 


APPENDIX  A 

AREA  AND  VOLUME  FOR  A CONE 


1 Calculation  for  area  in  radial  direction 


1.2  Arm 


Equation  to  calculate  average  cross  section  area  of  the  node  adjacent  to  the 
node  at  the  side  surface 


A r,  = - (2 it 

Ar2  = - (2n  f J H ,) 
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Equation  to  calculate  average  cross  section  area  of  the  node  at  the  surface 

Su  = JHJ  * 2 

2 Calculation  for  nodal  area  in  axial  direction 

2.1  Ayv 


Equation  to  calculate  average  cross  section  area  of  the  interior  nodes 
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Equation  to  calculate  average  cross  section  area  of  the  node  that  adjacent  to 
the  node  at  the  side  surface 


i McvlMipffjforp.9fe.8iym!; 


Equation  to  calculate  volume  of  the  interior  nodes  of  height  // 


APPENDIX  B 

GRID  GENERATION  IN  A FINITE  ELLIPTICAL  CYLINDER 


Figure  B.  1 An  ellipse 


An  ellipse  is  given  by  the  equation: 

ll  , = 


(B-l) 


where  ‘a'  is  one  half  the  length  of  the  major  axis,  'b'  is  one  half  the  length  of  the  minor 

A 2-dimensional  finite  difference  model  is  developed  for  a finite  elliptical  cylinder 
with  height  ‘h*.  The  finite  elliptical  cylinder  was  subdivided  in  to  series  of  concentric 
elliptical  ring  stacked  on  top  of  each  other  with  series  of  heat  flow  curves  drawn 
perpendicularly  with  the  ellipse  perimeters.  Due  to  symmetry,  only  1/4  of  the  finite 
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elliptical  cylinder  is  used  to  develop  the  energy  balance  equations.  Heat  flows  along  its 
height  and  along  the  heat  flows  curves  (Figure  3.13). 

The  input  parameters  for  generating  the  grid  for  a finite  elliptical  cylinder  arc 

• Dimensions  including  a,  b,  and  h 

• Number  of  the  concentric  ellipses,  number  of  the  heat  flow  curves,  and 
number  of  nodes  in  height 

Grid  generation  involves  determining  of  the  perimeter  of  an  ellipse  in  one 
quadrant  (P)  which  can  be  calculated  as 


and  E is  an  incomplete  elliptic  integral  of  the  second  kind  calculate  by  using  subroutine 
ELI2  (given  at  the  end  of  this  appendix). 

Once  the  quadrant  cord  is  found,  the  strategy  for  grid  generation  is  as  follows 
I Draw  a series  of  heat  flow  curves 

1 . 1 Locating  of  intersection  points  of  the  ellipse  and  the  heat  flow  curves 
The  quadrant  cord  is  divided  by  the  number  of  heat  flow  curves  to  obtain  equal 
arc-lengths.  By  using  bisection  method  and  subroutine  ELI2,  the  angle  from  the  center 


■2’ 


(B-2) 


where  ‘e’  is  an  eccentricity  of  the  ellipse 


(B-3) 


(6  )that  corresponds  to  the  equal  arc  length  can  be  found. 


Figure  B-2  Location  of  the  point  on  the  perimeter 


Form  Figure  B-2  the  radius  (r)  can  be  written  as: 


(B-4) 


The  coordinate  at  the  perimeter  is 


x - r cos  6 


(B-5) 

(B-6) 


1 .2  Draw  heat  flow  curves  from  the  center  to  the  intersection  points 
Once  the  intersection  points  on  the  perimeter  of  the  ellipse  are  found,  the  curves 
will  be  drawn  form  the  center  of  the  ellipse  toward  these  points  having  the  equation 

y-cx’  (B-7) 


These  curves  are  perpendicular  with  the  ellipse  perimeter  at  the  intersection 
points.  The  product  of  the  slope  of  the  ellipse  (SE)  and  the  slope  of  the  curve(SP)  equals 
- 1 . The  slope  of  the  ellipse  can  be  found  by  rearranging  ellipse  equation  (B- 1 ) as 
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(B-8) 


Tile  slope  ai  a point  on  the  perimeter  of  an  ellipse  (SE)  is  the  derivative  of  the 
equation  B-8. 


(B-9) 


The  slope  of  the  heat  line  curvc(SP)  is  the  derivative  of  the  equation  B-7 


(B-ll) 


SP  * SE  = -I  (B-12) 

The  power  of  the  heat  flow  curve  (n)  and  coefficient  (c)  could  be  solved  as 


(B-13) 


The  equation  for  desired  x by  the  Newton  method  is 


(B-I7) 


Once  the  x coordinates  are  found , the  y coordinates  could  be  obtained  from 
equation  B-8. 

Calculation  of  the  cross  section  nodal  area 


Figure  15-3  Illustrate  cross  section  area 

The  area  from  (a,o)  to  the  heat  flow  curve  at  of  intersection  point  (x,y)  as  shown  in  figure 
B-3  could  be  calculated  from 


(B-IS) 
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Sub  ELI2(R,  x,  CK,  a,  b) 


PURPOSE 

COMPUTES  THE  GENERALIZED  ELLIPTIC  INTEGRAL  OF  SECOND  KIND 
USAGE 

CALL  ELI2(R,X,CK,AB) 

DESCRIPTION  OF  PARAMETERS 
R -RESULT  VALUE 

X - UPPER  INTEGRATION  BOUND  (ARGUMENT  OF  ELLIPTIC 
INTEGRAL  OF  SECOND  KIND) 

CK  -COMPLEMENTARY  MODULUS 
A - CONSTANT  TERM  IN  NUMERATOR 
B -QUADRATIC TERM  IN  NUMERATOR 

REMARKS 

MODULUS  K = SQRT(I.-CK*CK). 

SPECIAL  CASES  OF  THE  GENERALIZED  ELLIPTIC  INTEGRAL  OF 
SECOND  KIND  ARE 

F(ATAN(X),K)  OBTAINED  WITH  A=  1 B=1 . 

This  is  what  we  want  E(ATAN(X),K)  OBTAINED  WITH  A=  I B=CK'CK 
k'2  + ckA2  = 1,  and  k~2  <=  1 

B(ATAN(X),K)  OBTAINED  WITH  A=l B=0. 

D(ATAN(X),K)  OBTAINED  WITH  A=0„  B=l 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
NONE 


METHOD 
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DEFINITION 

R=INTEGRAL((A+B*T*T)/(SQRT((1+T,T)*(I+{CK,T)**2))*(1+T*T)), 
SUMMED  OVER  T FROM  0 TO  X). 

EQUIVALENT  IS  THE  DEFINITION 
R=INTEGRAL((A+(B-A)*(SIN(T)),*2)/SQRT(1-(K*SIN(T)),*2), 
SUMMED  OVER  T FROM  0 TO  ATAN(X)). 

EVALUATION 

LANDENS  TRANSFORMATION  IS  USED  FOR  CALCULATION. 
REFERENCE 

R.  BULIRSCH,  NUMERICAL  CALCULATION  OF  ELLIPTIC  INTEGRALS 
ELLIPTIC  FUNCTIONS 

HANDBOOK  SERIES  OF  SPECIAL  FUNCTIONS 
NUMERISCHE  MATHEMATIK  VOL.  7,  1965,  PP.  78-90. 


Dim  C,  D,  AN,  AA,  ang,  PIM,  ISI,  ARI,  GEO,  SGEO.  AAR1,  AANG,  PIMA 

TEST  ARGUMENT 
If  x = 0 Then 
R = 0# 

Exit  Sub 
End  If 

TEST  MODULUS 
C = 0# 

D-0.5 

If  CK  = 0 Then 
R = Sqr(l#  + x * x) 

R = (a  - b)  * Abs(x)  / R + b * Log(Abs(x)  + R) 


TEST  SIGN  OF  ARGUMENT 


R = b 

ang  = Abs(  1 # / x) 


GEO  = Abs(CK) 

' LANDEN  TRANSFORMATION 
Landen: 

R = AA  * GEO  + R 
SGEO  = ARI  • GEO 
AA  = AN 
AARI  = ARI 

' ARITHMETIC  MEAN 
ARI  = GEO  + ARI 
1 SUM  OF  SINE  VALUES 
AN  = (R  / ARI  + AA)  * 0.5 
AANG  = Abs(ang) 


PIMA  = PIM 

If  ang  = 0 Then  ang  = -0.0000000 1 
If  ang  <=  0 Then 


PIM  = PIM  + 3.1415927 
ISI  = ISI  + I 
End  If 

AANG  = AR1  * AR1  + ang  * ang 
p = D/Sqr(AANG) 

If  ISI  >=  4 Then  ISI  = ISI  - 4 
If  ISI  >=  2 Then  p = -p 
C-C+p 

D = D * (AARI  - GEO)  * 0.5  / ARI 

If  Abs(AARI  - GEO)  > 0.0001  • AARI  Then 
SGEO  = Sqr(SGEO) 

• GEOMETRIC  MEAN 
GEO  = SGEO  + SGEO 
PIM  = PIM  + PIMA 
ISI  = ISI  + ISI 
GoTo  Landen 
End  If 

1 ACCURACY  WAS  SUFFICIENT 
R = (Aui(ARI  / ang)  + PIM)  • AN  / ARI 
C = C + D * ang  / AANG 
R=R  + C*(a-b) 

If  x < 0 Then  R = -R 
End  Sub 


APPENDIX  C 

EXPERIMENTAL  AND  SIMULATED  TEMPERATURE  IN  AN  ACRYLIC  MATERIAL; 
FINITE  ELLIPTICAL  CYLINDERS  AND  CONES 


B\p  I l \p  2 Sim  I 


79.03  37.21  26.45  35.94  26.78 


79,61  53.58 


using  h = 1300  W/  in2  K.  K = 0. 195  W/m  K,  Alpha  =1.1 


Table  C-3  Experimental  and  simulated  data  for  the  finite  elliptical  cylinder 
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Exp  I E\p  2 I '!'  ■ I \ [ ■ i 


Sim  2 Sim  3 Sim  4 


73  38.98  47.2 


50.87  59.37  71 
54,49  62.54  72 


1.62  9.21  15 

75.56  12.71  15 


74  54.52  63 


53.29  62.37 

43.19  58.48 
39.11  55.85 

57  32.3  50.3  42 


Tabic  C-4- 
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64.65  70.74  75.01 


75.9  71.83  61.43 


95  69.41  72 


24  76.35  75 


97  75.58  72 


number  of  ellipse  =16,  number  of  heal  flow  curve  = 1 6,  and  node  in  height 
This  is  llic  interporated  data  lo  show  temporature  for  every  30  second 


for  the  cone  # 1 , Experiment  I 


42.2  29.43  28,82 


65.08 

66.59  73 


93  67.73  73 


63  71.27  75,91 


72.59 


Tabic  C-5 
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I l \p  1 Exp  2 


85.72  73.83  73, 

75.25  75 


77.82  24.86  42.76  43 


'6.75  72.7  56 

83  75.21  68.01 


95  47.29  47 


32.77  43.61  43.22 


Tabic  C-6  Experimental  and  simulated  data  for  the  cone  4 1 . Experiment  2 


Water  Exp  I Exp  2 Exp  3 Exp  4 Sim  I Sim  2 Sim  3 Sim  4 
25.34  25.16  25.13  25.08  24.7  25.02  25.02  25.02  25.02 

25.26  31.11  25.11  25.02  25.23  31.29 

16.15  38.53 


50.59  48.71  55 


23  55.45  56 


52.7  57 

54.53  59 


Table  C-7  Experiment  and  simulation  data  for  die  cone  ~ 2,  Experiment  I 


35  67.84  75.99 


n!  K,  K = 0. 19  W/m  K,  Alpha  = 1 .09  xIO'7  m2/sc 


Table  C-8 


Time  Tenipcralurc  (°C) 


23,93  36.11 
24.7  38.18 
25,55  40.1 


55.81  50.22 

56.52  51.09 


31.79  48.27 


42.36  34.09 


35.25  51.58 


712  6975  72  X6 


n!K,  K = 0. 1 9 W/m K.  Alpha  = 1.09  x 1 O'1  it 


This  is  the  interporated  data  to  show  temperature  for  even'  minutes 


Table  C-9  Experimental  and  simulationcd  data  for  the  cone  H 3,  Experiment  I 


Exp  1 Exp  2 Exp  3 


29  26.27  25 


25.59  25.5 


05  39.72 

42.6  29 


57.52  45.2 
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(min)  Water  Exp  I Exp  2 Exp  3 Sim  1 Sim  2 Sim  3 


APPENDIX  D 

SIMULATED  HEATING,  COOLING  AND  TOTAL  PROCESSING  TIME 
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APPENDIX  E 

SIMULATED  QUALITY  RETENTION 


Tabic  I:'-  I Simulated  quality  retention  of  foods  with  thermal  diftusivity  1 .2  x 10'7  m:  /s 
after  thermal  processing  in  the  conically  shaped  container  #1,  h = 300  w/m‘K 
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Table  E-  2 Simulated  quality  retention  of  foods  with  thermal  diffiisivity  1.4x10*  m*  /s 
after  thermal  processing  in  the  conically  shaped  container  41,  h = 300  w/m2K 
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Table  E-  3 - continued. 


Table  H-  4 —continued. 


able  E- 
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Table  E- 1 


Tabic  E-7--cominucd. 


Table  E-  8 Simulated  quality  retention  of  foods  with  thermal  diffusivity  I A x I O'7  in:  s 
after  thermal  processing  in  the  conically  shaped  container  32.  h = 300  w/nrK 


Tabic  E-8--contimicd. 
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TableE-9  Simulated  quality  retention  of  foods  with  thermal  diflusivity  1.6x  10'  ill  /s 
after  thermal  processing  in  the  conically  shaped  container  tt2,  h = 300  w/mJK 


D (min)  Temp(°C) 


Table  E- 9~o 


Table  E-10-c 
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Table  E- 
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Table  E-  13  Simulated  quality  retention  of  foods  with  thermal  diffusivity  1.2  x 10'  m /s 


after  thermal  processing  in  the  conically  shaped  container  #3,  It  = 300  w/m3K 

% Retention 


Table  E- 


14  Simulated  quality  retention  of  foods  with  thermal  diffusivity  1.4  x 10*'  m"  Is 
after  thermal  processing  in  the  conically  shaped  container  #3,  h = 300  w/m’K 


D (min)  Tcmp(°C) 


Table  E-  14-continued. 


Table  E- 


% Retention 

D (min)  Temp(°C)  Z = l(l"C  Z-2tl°C  Z = 30°C  Z = 4(I°C  Z«50"C  Z-60°C 


Table  E- 19  Simulated  quality  retention  offbodswilh  thermal  diflusivity  l.2x  10'7  /s 
after  thermal  processing  in  the  clliptically  shaped  container  #1,  h = 300  w/m!K 


Table  E -19-continued 


Table  E-20  Simulated  quality  retention  of  foods  with  (hernial  diffiisivity  1.4*  lO’m’/s 
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: E-  21 
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Table  E-22  Simulated  quality  retention  of  foods  with  thermal  diftusivity  1 2 x I0*7  nr/s 
after  thermal  processing  in  the  elliptically  shaped  container  Hi  using  Infinite  h 

% Retention 


D (ntin)  Templ°C) 


is?  ’■asp  “ar 

!KS  E S SS  5S 


ss  s 
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after  (hernial  processing  in  (he  elliptically  shaped  container  #1  using  Infinite  h 
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Table  E- 23  -continued 
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Tabic  E-  24  Simulated  quality  retention  of  foods  with  thermal  difiusivity  l,6x!0'7m2/s 


' t™pm:  Z»ZC  \-'4".;c 

is  ss  s s as  s as 

S X S jffi  S:S  S Ej 

is  g § g Is  If  g 

S g Ijg  S3  S S|  X 

E IE  II  g g g Is 

is  E s g E Si  g 

;g  s ss  s s s s 

E IS  If  g If  5 g 

ii|  g g If  g §g  g 

E g g g if  g g 

ill  g g g g g | 

E g II  II  IS  II  g 


Table  E-  25  Simulated  quality  retention  of  foods  with  thermal  diflusivity  12  x 10  7 is 
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Table  E- 26  Simulated  quality  retention  of  foods  with  thermal  diffusivity  1.4x  ! 0*7  m'  /s 
after  thermal  processing  in  the  elliptically  shaped  container  #2,  h = 300  w/m!K 
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